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PREFACE 


I don’t know whether many people ever bother to read prefaces 
but I feel that anyone who wishes to find what a book is about should 
at least be offered the opportunity at as early a stage as possible. 

First I must make it clear that, although I hope there is something 
to be learned from this one, it is in no sense a ‘teaching’ book: 
its intended purpose is to interest potential mathematicians and to 
introduce them to a branch of arithmetic which may transform their 
schoolday impressions of the subject. 

Curiosity is one of the few virtues common to all scientists but 
unfortunately today its satisfaction is mostly dependent upon expen- 
sive equipment and closely knit team work. The research laboratory 
of a Theory of Numbers enthusiast can, at a pinch, be equipped 
with nothing more than pencil and paper. And to anyone who 
hesitates to set up such a laboratory and embark upon his own 
research because of the view that it has all been done before and 
there is nothing left to discover I would say with emphasis that 
this is not true. In mathematics there is always another field beyond 
the hedge if one can only find the gap to crawl through. 

The main themes of this book are ‘Factorisation’ and ‘Prime 
Numbers’ but there are some diversions into such topics as recurring 
decimals, additive series, a brief note on Diophantine equations 
and a few others. They have all been chosen because they are among 
those branches of Number Theory which I have found intensely 
interesting—and I hope some of this enthusiasm will rub off during 
the reading—and also because they show how unexpectedly some 
divergent lines of enquiry latch in together. 

I have included a chapter (Finite Arithmetic) on elementary 
congruence technique since some knowledge of this brilliant inven- 
tion is essential to any serious study of large numbers. For the same 
reason there is a brief note on the ‘Converse of Fermat’s Theorem’ 
and its relation to composite numbers. 
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A number of tables of values and functions have been included in 
the text and the Appendix; these have the double service of providing 
examples for the former and at the same time saving the potential 
research worker a certain amount of muscular mathematics. 

Most readers will, I am sure, find something about the properties 
of numbers which is new to them in the following pages. (How 
many professionals, for instance, know of a direct linkage between 
1/79—having a recurring period of thirteen digits—and +/2 which 
is irrational?) 

However it is the amateur mathematician I have had most in 
mind in compiling the following notes, and in defining my concep- 
tion of the term ‘amateur’ I cannot do better than quote J. E. 
Littlewood (A Mathematician’s Miscellany, Methuen & Co. Ltd.). 

Ἵ constantly meet people who are doubtful, generally without due 
reason, about their potential capacity. . . If your education just 
included, or just stopped short of including, ‘a little calculus’, 
you are fairly high in the amateur class.’ 

Coming from such a master this is pretty heartening and I hope 
his words will encourage those who find enjoyment and relaxation 
in solving puzzles to apply their ingenuity to experimenting with 
numbers. 

Finally I would like to stress that in studies of this nature it is 
particularly important that the reader takes more than a passive 
interest in the text and I have therefore ended each chapter with a 
few relevant exercises. (The answers will be found at the end of the 
book.) 


I am greatly indebted to Mr. Nicholas E. Scripture for some 
valuable advice and a critical scrutiny of the manuscript. 


COUNTING 


(1) In this book we shall be concerned almost entirely with the 
integers. Those all too familiar numbers 1, 2, 3,. . . which go on 
interminably, each one formed by adding ‘one’ to the last. 

When man first began to feel the need for more specific terms 
than his current equivalent sounds for ‘a’, ‘few’ and ‘many’ he 
began to count. And counting is the basis of all arithmetic for it 
involves the process of addition. 

One can only speculate on the evolutionary processes that led up 
to the system of counting by the addition of units as we know it, 
or on the many experiments which must have been tried before we 
settled down to counting generally in bundles of ten units. It has 
been plausibly put forward that the decimal scale originated in our 
having ten fingers on which to count but the weights and measures 
and the monetary systems of the world give this view little support. 

Be that as it may we now have a stable and familiar arithmetic 
which uses the ten symbols, 0, 1, 2,. . . 9, and expresses numbers 
greater than nine in a form of shorthand based on the powers of ten. 
Thus the number ‘two thousand five hundred and seventy nine’ 
contains 9 x 100 - 7 x 10!+5 x 10?+2 x 10° units and is 
conventionally written 2579. 

(In the power index notation n° τὸ 1, and πὶ =n.) If 2579 is 
successively divided by ten the remainders 9, then 7, 5, and 2 
are left, these being its own digits reading from right to left. 

All very elementary of course, but only because we are so accus- 
tomed to this system that there is no longer any need to stop and 
think how the numbers we use are constructed. Their expression as 
multiples of powers of ten is, as we shall see, well chosen but it is 
none the less a matter of pure convention and it should not be 
forgotten that the same numbers could just as well be expressed as 
multiples of powers of any base we care to choose. 

In order to convert a number in the decimal scale to that of any 
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other base it is only necessary to divide successively by the number of 
the new base writing down the remainders in the reverse order to 
that in which they appear. Converting the number 2579 into the 
scale of seven for instance we have: 


Remainder 
7)2579 3 
7) 368 4 
ἜΣ ἊΣ 3 
ἥν. 0 
Ι 1 


the new expression now being written 10343, indicating that it is 
equal to: 


3x 72 =e3 x l= 3 
4x7i=4x T= 28 
3x 7=3x 49= 147 
0x 72= 0 
1 x 74*=1 x 2401 = 2401 

2579 


It will be seen that in dividing by seven there can be no remainder 
greater than six and so the only digits that appear in this scale are 
0, 1,2,. . . 6. Similarly in the scale of two we require only the two 
digits Ὁ and 1. 

For example: 


», ὦ “ ocoooorcoor = 
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the decimal number 2579 now appearing as 101000010011 which 
states explicitly that it is equal to 
29 + 21 + 24 4 29 4 QU 
In just the same way the number could be given in powers of any 


number greater than ten, say for example twenty-three. (The divi- 
sions by 23 are given here in condensed form.) 


23)2579 3 
23) 112 20 
4 4 


In writing this down we have to remember that the ‘20’ must 
be regarded here as a single digit; this can be made clear by enclosing 
it in brackets and the new number now appears as 4(20)3 indicating 
that 2579 is also equal to 3 + (20 x 23) + (4 x 23%). 

It will be convenient to denote any arrangement of digits as Ng, 
in which s represents the scale of notation employed: the observa- 
tions that have been made can then be written concisely as follows, 


(N10) 2579 = (Nz) 10343 
= (N2) 101000010011 
= (Nes) 4(20)3 


It is now seen that whilst the choice is quite arbitrary the ‘powers 
of ten’ notation presents a good all-round middle course for express- 
ing numbers requiring more than one digit; the use of fewer symbols 
calls for longer strings of digits to represent a given number while 
more symbols would have meant more and longer multiplication 
tables to memorise. 


(2) There are occasions, however, when scales other than the decimal 
can be used with advantage. One of these will have become familiar 
to readers through the widespread development of electronic 
digital computers. Here the binary system which employs the two 
digits 0 and 1, is readily correlated with the two (on-off) positions 
of an electrical switch whilst an increase in the number of digits is of 
small consequence. 

Again, because the binary system’s largest digit is 1, it follows 
that any number whatever can be expressed as the sum of single 
powers of 2. Thus it is a simple matter to multiply together any 
two numbers without going outside the x2 multiplication table. 
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This is done quite easily by successively multiplying one of the num- 
bers and dividing the other by two and then deleting those numbers 
in the increasing column which correspond with the even numbers 


in the other. 


In effect we are thus finding by successive division those powers 


of two which are needed to represent one number, multiplying the 
other by 2, 22, 23, etc., and then rejecting the products which are 


not required. 
For example: To find 37 x 127(= 4699). 
Divide by 2 Multiply by 2 

37 | 127 
18 O——_ ir 

9 | 508 

4. QO— HC 

2 0---------------2032- 

I | 4064 


4699 = Answer. 


—_- - 


The left hand column will be familiar; in it we have converted 
(Nio) 37 to (Ne) 100101 showing that 


37 = 294 224 28 = 14+ 44 32. 


The right hand column is made up of 127 multiplied in turn by 
5, 5, Rep .ὄὯ5, and from it we have deleted the products of 
2, 23, and 24, the remaining figures then being added. In practice 
there is no need to write down the first column remainders since the 
zeros occur alongside even numbers and these can consequently 
be used to locate the deletion lines. 

The process is much simpler to carry out than to describe and it is 
rather surprising that it is not taught as a rule of thumb method in 
those schools whose children have difficulty with their multiplication 


tables. 


(3) Apart from these instances one might readily conclude that there 
is little to be gained by any further study of scale-changing and in- 
deed this is about as far as most textbooks take the subject. Never- 
theless it is never safe to assume that any branch of mathematics 
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has arrived at a dead-end and it will be seen that further develop- 
ments are not only possible but of practical interest. But first it 
will be necessary to make a digression over well trodden ground. 
There are two well known and simple tests for determining whether 
a given number is divisible by nine or eleven. For the former the 
digits of the number are added together to form another number 
and if this contains more than one digit the process is continued 
until only one remains. If this digit is 9 then the original number is 
exactly divisible by 9; any other resultant digit will be the remainder 
found on dividing the original number by 9. Thus the number 
23456723 when divided by nine will leave the remainder 5, since 
2+34+44+54+64+74+2+43 = 32, and 342 = 5, The rule 
commonly known as ‘casting out the nines’, can also be applied to 
— 4 three, the remainder in the above case then being 
Somewhat similarly a number is exactly divisible by eleven if the 
pees gh = sums of its alternate digits is 0 or a multiple 
of eleven. For example 73645 is divisible by eleven si 
a eg y eleven since (5 + 6 + 7) 
Now the numbers nine and eleven are respectively 10 — 1, and 
10 + 1 and it would be reasonable to expect that the same rules 
apply to the divisors s — 1 and s + 1 when dealing with numbers 
expressed in the scale of s. Suppose for example we take the number 
4571 (= 7 x 653) and transpose it into the scale of 8, thus: 


ον ὅν... 3 
571 3 

71 7 

8 0 

1 1 


a (N10) 4571 = (Ng) 10733. 
nd since in the scale of 8 the ‘nines’ rule applies to division b 

7 
we have 3+3+7+1=14=2x/7, and th s 
divisor of 4571. ) a ok: 


(We could have continued, 14 = (Ng) 16. And 6 + 1 = 7) 


Or again, converting to the scal i 
, e of six and using the ‘ : 
rule we have: : ὌΠ 


(N10) 4571 = (Ne) 33055 


EEE 
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and as 5+0+3=5+3 it is shown that 4571 is divisible 
by (6+ 1) =7. 


The principle is in fact quite general and can be proved without 
difficulty. To take, for simplicity’s sake, a specific example, in the 
decimal scale every number can be expressed as follows: 


a+ 10b + 100c + 1000d+...=a+54+9b+c¢4+99¢+d+ 
999d+...=a+b+ec+d+...+ 9m (a multiple of 9) 
That is, the sum of the digits a, b, c, etc. will be the remainder after 
dividing by nine. 


The ‘eleven’ rule and a generalisation into any scale of notation 
can be deduced in an exactly similar manner. 


(4) We can now begin to glimpse another practical use for the 
process of scale changing. In many mathematical operations it is 
often desirable to know whether a given number is prime, or if not 
what are its factors. Provided the number is not too large the straight- 
forward approach is to make tentative divisions by the successive 
primes 7, 11, 13, 17, etc., until a factor is found or the square root 
of the number is reached. (Obviously if there is no factor less than 
the square root there cannot be one which is greater and the number 
must be prime.) From what we have seen above it will be clear that 
when working through the primes in this way there will occur many 
cases when two tests can be carried out by making a single scale 
change. Thus whenever two primes differ by two, that is when they 
are of the form ἢ — 1 and n + 1 they can be tested simultaneously 
by changing the given number into the scale of m and applying 
the ‘nine’ and ‘eleven’ rules to the new digits. 

Let us suppose for example that in attempting to factorise the 
number 18383 all the primes up to and including 23 have been 
tested without a factor appearing. Now instead of continuing with 
trial divisions first by 29 and then by 31 we convert the number into 
the scale of 30, thus: 


30 [[ ~=—s: 18383 23 
612 12 
20 20 
Applying the two tests to the ‘digits’ of the new scale we have, 
(a) 23 + 12 + 20 = 55 
(0) 23 + 20 — 12 = 31 
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Since 55 is not a multiple of 29 then this is not a factor but the 
‘eleven’ rule shows that 31 is. Dividing out we find that 18383 = 
31 x 593 and as 31? is greater than 593 the latter is also prime. 


(5) Now whilst it is possible to test both the prime divisors 41 and 
43 in one operation, the conversion into the scale of 42 is ‘awkward’ 


- and it is doubtful whether much time would be saved. The method 


is clearly most effective when, as in the last example the divisors to 
be tested lie on either side of a multiple of ten (i.e. 10” + 1). 

But, just as ‘casting out the nines’ can be used to detect multiples 
of three—a factor of nine—the above method is equally applicable 
to primes which are factors of numbers of the forms 10” — 1 and 
10n + 1. For instance the scale of fifty provides an immediate 
test for the divisors 7 and 17 because 49 = 75, and 51 = 17 x 3. 

Taking for example the number 12733 (which equals 7 x 17 x 
107) and converting to the scale of fifty, we have: 


50 12733 33 
254 4 
5 5 


and we see that 


33 - 4 -- 5 τὸ 42 =amultiple of 7 
33 + 5 — 4 = 34 = a multiple of 17 


This property leads at once to an extremely simple test for prime, 
or prime factors of, numbers of the form 10" + 1. In these cases no 
actual division is required to make the scale change as, for example 
1317459 = (Niooo), 1(317)(459). Now 999 is a multiple of 37, 
and 1001 = 7 x 11 x 13 and 


since 459 + 317 + 1 = 777 = 37 x 21 
and 459 + 1—317 = 143 =11 x 13 


it is seen at once that 1317459 is exactly divisible by 11,13 and 37, 
but not by 7. 

In general a given number can be tested for divisibility by any of 
the factors of 10" + 1 by separating it into sets of nm digits from 
the right and then applying the ‘nine’ and ‘eleven’ rules to these sets. 

Incidentally it will now be seen that the basic test for, say, divisibi- 
lity by eleven is by no means the only one. In fact, since the successive 
numbers 10 + 1, 100 — 1, 1000 + 1, 10000 — 1, and so on are 
all multiples of eleven it is a simple matter to devise a test tailored 
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to suit the magnitude of any given number. Briefly, the number is 
divided into large ‘sets’ and then, according to the number of digits 
per set the ‘nine’ or ‘eleven’ rules are applied until only a two digit 
number remains. If this is divisible by 11 then so is the original 
number. 

Thus if we wish to find the remainder when the number 7970684847 
is divided by 11 we proceed first with subtraction because the 
number can be separated into sets of five digits: 


84847 
79706 
5141 
and then by addition of pairs of digits: 
51 
41 


2 = (11 x 8) - 4. 
And therefore the remainder is 4. 


(6) It will be convenient to draw up a table of ‘easy’ scale-changes 
showing the prime numbers to which the (s—1) and (s+ 1) 
rules can be applied as tests of divisibility. By ‘easy’ I mean those 
scales which do not in effect require successive division by any 
number greater than eleven. 

After what has been said the following Table 1 should be self 
explanatory; the numbers in brackets are all primes whose divisi- 
bility can be tested by changing the given number into the scale ‘s’, 
and it will be seen that all the primes less than 100 are covered with 
the exception of 73, 83, and 97. And with some of them a single 
scale-change may eliminate as many as four possible prime factors. 

Provided we remember which column we are working in it is 
quite practicable to change horses in mid-stream. Suppose, for 
instance, it is required to find the remainder, if any, on dividing 
123456 by 17, then we have, 


800 | 123456 256 
154 154 


410 50 | 410 10 
8 8 


2 = remainder 


The exceptions mentioned, namely 73, 83 and 97, could be 
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Table 1 


Scale(s) 


299 (13, 23) 


799 (17, 47) 
899 (29, 31) 
999 (37) (7, 11, 13) 


dealt with in the scales of 220,250 and 290 respectively but these 
do not present such simple divisors as the above. 


EXERCISES 

1. Express the integer 11111 in the scale of 2. 

2. Express the integer 11111 in the scale of 16. 

3. Express the integer 11111 in the scale of 32. 

4. What observations can you find to make on the answers to 
exercises 1, 2, 3,? 

5. Multiply 123 by 456 by the process shown in the text and check 
by reversing the procedure. 


6. Without using ordinary long division find what remainder, if 
any, is left after dividing N = 455331 by 707. 
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(1) In the last chapter we saw how it was possible to tell whether a 
given number was a multiple of certain primes by dividing it, in 
effect, by a smaller number. The method was developed by 
generalising two well-known arithmetical tricks. Let us look at 
these two again from another angle. 

Numbers in the decimal scale which are exactly divisible by nine 
are recognisable at once, as shown in /, 3), from the fact that the 
sum of their digits is always nine or some multiple of nine, thus: 

2x9=18 and 8+1=9 
3x9=27 and 7+2=9 
4x9=36 and 6+3=9 
29773 Χ 9 = 267957 «=and = 74+54+9+7+6+2= 36 
and so on. 

The above test deals with the number 10¢— 1, where ¢ = 1. 
It would be natural to expect that similar rules could be devised for 
those numbers where ¢ takes the values 2, 3, 4,. .., and it is 
readily seen that this is indeed the case. 

Taking, for instance, some multiples of 19, where ¢ = 2, and 
operating upon them by first multiplying the ‘units’ digit by two 
and then adding the result to the ‘tens’ digit it is found that the new 
number is then either 19 ora smaller multiple of 19. In the latter case 
the operation is repeated as in the method of ‘casting out the nines’. 

38 57 76 247 84493 (= 19 x 4447) 
16 ie 12 14 bails. 
19 19 19 38 8455 


16 10 
19 855 
10 
95 
10 
19 
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Similarly when ὦ = 3, and it is desired to test for multiples of 29, 
the multiplier now becomes 3 and we have: 


58 87 3393 12673 (= 29 x 437) 


2 42] gees 
oS Ὁ 348 1276 
24 18 
58 145 
24 15 
29 29 


The principle is clearly general and can perhaps best be explained 
by re-examining the procedure when ¢=2. What we did, in . 
effect, was to multiply the units digit by 19, or (20 — 1), add, 
and then divide by 10. 


For example 38 
152 = 8 x 19 
19 
Now let us suppose that N = 19m = a + 10b. Then, adding 19a, 
we have 19a +a+ 10} = 102a+ δ). Since 10 is not a divisor 
of 19 then 2a + ὃ must be another multiple of 19 less than 19m 
and if the process is continued we must eventually arrive at the 
smallest possible multiple, namely 19 itself. 
Looking again at the examples given above it will be seen that 
when WN is exactly divisible by 19 the process discloses the quotient 
m. For instance, in the case of N = 247 we have: 


247 247 = 19m 
14 | 133 =19x7 
38 is ‘shorthand’ for 380 
16 1520 = 19 x 80 
19 1900 = 19 x 100 
Therefore 19m + (19 x 7) + (19 x 80) = 19 x 100 
or m = 100 — 87 = 13 


In the algorithm the underlined digits, 87, provide us with m 
after subtraction from 100. 
Similarly when N = 84493 the terminal digits show that since 


10000 — 5553 = 4447, then 84493 = 19 x 4447. 


la at I AO CY OE AS Re σὸν 
---͵ .Ρ- 


~ ee EE .τπὰνἱἷς .-π' τὶ 


22 Casting out the Primes 


Or again, where we were dealing with multiples of 29, in the case 
of N = 12673, the terminal digits are 563. 
Then 1000 — 563 = 437 and hence 12673 = 29 x 437. 


(2) It is important at this stage to note that nothing has been said 
to imply that this method of division only has its uses when 10t — 1 
is a prime. In fact.it can be applied with the same advantage to any 
of the prime factors of 10f—1, as the following examples will 
show. 

Since 13 x 3 = 39, wecan use 4 as the multiplier to show that 
4771 is a multiple of 13 thus: 


4771 
4 
481 
4 
52 
8 
13 

And, multiplying the terminal digits 211 by 3 and subtracting 
from 1000 we have the quotient 367. 

Similarly, noting that 13 x 23 = 299, it will be clear that, using 
the multiplier 30 as above, it is possible to test both these factors 
simultaneously. 

Given N = 12029 (= 23 x 523) then, 

12029 
270 
1472 
60 
207 
210 
23 


Thus 23 is a factor of 12029 but 13 is not. The quotient can be 
found as in the previous example. 


(3) It is more than likely that the reader is already anticipating the 
next step in developing these principles. If this is so it might be a 
good thing for him to close the book at this point—temporarily 1 
hope—and conduct some experiments of his own. 
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Whether this suggestion is acted upon or not our experience 
of the effects of scale-changing make it imperative to examine the 
treatment of those divisors which have, or can be put into, the form 
of 10¢ + 1. As one might expect the method of test depends again 
upon the value of ¢, and it is quickly seen that by multiplying by 
t, etc., as in 2.1, but in this case subtracting, any multiple of 
10: + 1 finally yields the digit 0. 

Thus when t = 2 we have: 


42 105 30639 (= 21 x 1459) 


4 10 18 
2 oe 
10 
294 
8 
21 
2 
0 


What is more, as the last example shows, the terminal digits 
1 4 5 9 indicate the quotient directly. 

As before the test is equally applicable to the factors of 10¢ + 1 
when this number is composite. In such cases, however, the method 
may produce a small multiple of the prime under test instead of 
zero, and the quotient is then rather more elusive. In order to keep 
this chapter within reasonable bounds I shall give only one example 
of this particular condition and the reader will be left to discover a 
general process for finding the quotient. 


Let N = 321433 (= 47 x 6839). Now 47 x 3= 141, and we 
shall use the ‘multiplier’ 14, thus: 


321433 
42 


32101 
14 


3196 
84 
235 
70 
—47 


...-.-...»»» 
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This example has been given not only to illustrate the above state- 
ment but also to show that in many cases this method of testing 
possible divisors of N has a distinct advantage over the scale- 
changing technique. As we saw in Chap. 1 the latter is at its best when 
it can be applied to primes related to what I have called ‘easy’ 
divisors. On the other hand, in this system however large we have to 
make t¢ it has only to be multiplied by the terminal digits, none of 
which of course can be greater than 9. 


(4) As we shall see later, when it comes to the factorisation of 
large numbers it is important to eliminate as many as possible of 
the small prime factors before starting on more sophisticated 
techniques. In the ordinary way this is a laborious process and any 
method of shortening the work is to be welcomed. The algorithm 
we have been experimenting with provides such a method and it 
will now be convenient to gather together our findings into a form 
available for easy reference. 

To test whether a given number is exactly divisible by a prime p, 
choose a suitable multiple of p such that 


mp=10t-—1 or 10t+1 


CASEA Where mp = 10t — 1. 

The number to be tested (NV) will, in the decimal scale be of the 
form N= 10b+ a. From this derive the number N, = ὃ + fa. 
Νι then equals 10d + c, and again from this derive No = ὦ + te. 
Continue in this manner and if eventually a number ΝᾺ is reached 
which equals either p or a small multiple of p then p is an exact 
divisor of N. 


CASEB Where mp = 10t + 1. 

Again given that N= 105+ a, the numbers Ni = ὃ — fa, 
Nze=d-—tc, and so on are formed successively as above. If 
eventually it is found that Nz = 0, or p, or a small multiple of p, 
then p is an exact divisor of N. 

This sounds, and probably is, an unnecessarily formal way of 
describing a process which is extremely easy to operate in practice. 
This is a difficulty which is bound to arise when one employs methods 
which do not appear in ordinary arithmetic and is the principle 
reason why so many examples have been presented. 

Now whilst m can be chosen so that pm is equal to either of the 
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numbers 10¢ + 1 it will be obvious that for the purpose we have in 
mind one of these will be preferable to the other. For instance, if 
we wish to test whether 23 is a divisor then we have 23 x 3 = 69, 
and t = 7. To finda multiple of 23 of the form 10¢ + 1 we would 
have to take 23 x 7 = 161, 1 then becoming 16. Neither of these 
values presents any difficulty but I think most people would choose 7 
as the multiplier and the ‘addition’ process for the test. 

In the following table the primes below one hundred are listed 
together with their appropriate testing multipliers “(+) and «(—). 


Table 2 


(5) The methods so far described have been developed from ele- 
mentary properties of numbers and are in fact nothing more than 
extensions of the well-known rules for ‘casting out the nines’ and 
detecting multiples of eleven. Constructing algorithms of this sort is 
always an interesting—and often rewarding—form of mathematical 
recreation and it is surprisingly easy to devise divisibility tests from 
purely empirical observations. 

For instance the fact that 17 x 6 = 102, immediately suggests 
a simple method of dividing by 17. One has only to double the 


ial = = 
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left-hand digit, subtract from the next pair of digits and continue 
in this manner, thus: 
N = 13579 (= (17 x 798) + 13) 
13579 
2 
379 
6 
319 
6 
13 = the remainder. 
Multiplying the underlined digits 133 by 6, (102/17) gives us 798 
which is the quotient. 
We could, of course, have paired the digits off from the right and 
proceeded : 


ὃ SE Ἢ 
2 
Ὦ Ἢ 

66 
13 


This system is merely long division in a slightly different dialect 
but it is clearly preferable when dealing with large numbers and 
divisors which are closely adjacent to powers of 10. Thus one of our 
earlier ‘awkward’ divisors, 97, presents no difficulty when it is 
treated as (100 — 3). Here we use the multiplier 3 and add, as in 
the example, 


N = 1234567 (= (97 x 12727) + 48) 
1234567 


48 = the remainder 
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And again the underlined digits provide the quotient. Obviously 
it is a very simple matter to divide by any of the factors of numbers 
of the forms 9999... 99, and 1000. . . 01. 

I have yet to meet a “Theory of Numbers’ addict who can resist 
a problem and to round off this chapter I present the following little 
curiosity. Again use is made of the “ἢ function, this time in the form 
of the rising powers 10, 11, 7,. . .2, (note that = 1, 11} = 2), 
and the two worked instances should make the principle plain without 
the need of verbal explanation. 


(A) N = 31711 = 19 x 1669. = (¢ = 2) 
3x l= 3 
1x 2= 2 
7x 4=28 
1x 8= 8 
1 x 16 = 16 
57 
5x l= 5 
7x 2=14 
19 
(B) N = 34561 = 17 x 2033. (( -Ξ 5) 
2% im 3 
4x S= 20 
5x 25 = 125 
6 x 125 = 750 
1 x 625 = 625 


Now applying the ‘eleven’ rule we have, 
(20 + 750) — (3 + 125 + 625) = 770 — 753 = 17. 
The point of interest is that the increasing powers of ¢ are applied 
in the opposite direction (that is, from left to right) to what one 
might expect. 
EXERCISES 

1. If 43. = 21010101, find m by two different calculations, but 

without using ‘long division’. 
2. What remainder is left after dividing 110999999 by 997? 
3. Devise a simple process for dividing by 167. . 


3 RECURRING DECIMALS 


(1) If we take a large multiple of ten, say 10% where n> 20, 
and divide it in the normal manner by 19 we produce the sequence 
of digits 5263157894736842105263 ... which repeat themselves in 
bunches of eighteen digits. Dividing this number by the appropriate 
10” we are then able to say that 1/19 = -05263157894736842i, 
the suffixed periods indicating that this decimal part is repeated over 
and over again. 

There is a simple but important distinction to be observed here. 
The decimal -9999..., in which the 9’s are continued inde- 
finitely, and usually expressed as -9, is indistinguishable from unity 
and the above decimal -052... 842i is an exact divisor of “ὁ. 
On the other hand the number 526. . .421 is an exact divisor of 
999,999,999,999,999,999. Similarly 142857 is a divisor of 999,999 
and hence -142857 is a divisor of -9: it is in fact equal to 1/7. 

Reverting again to the fraction 1/19, we saw in the last chapter 
that division by 19 can be done, in effect, by multiplying succes- 
sively by 2 starting from the right. It will be seen at once that, 
starting with the digit 1, the above decimal is an excellent example 
of this process. 
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Thus, following this procedure we have: 
168421 


4210526315789473684210 


which, continued indefinitely produces the decimal recurrent cycle 
of the fraction 1/19. Another way of expressing the above summa- 
tion is 20° + 20! + 202 + 203+. .., which of course can be 
summed as follows, 


. « «368421 


The above sequence of digits can be arrived at in a number of 
ways. Two of these will be illustrated next in order to indicate the 
diversity of approach which this subject provides. 

The following tabulation will probably demonstrate one method 
more clearly than any laboured description of the process. 

Three columns, A, B, and C are set up in which B, is the 
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integral part of 4,2 and Ca is the remainder. The A, entries, 
starting with 1, are found from the relationship 


An +1 = 10Cn + Bn 


A B Cc 
2into 1 = 0 remainder 1 
10 

5 

12 


\o 
Ore NROTHA wWwAIPWOTNUK WAN WY 
mooocrrK OK ORR eS KF OOF SO 


Column B then provides the required sequence of digits. 

The next example is again best expressed in columnar form. This 
time the successive A entries are multiplied by 10 and we have the 
relationship An +1 = 10C, — B, where B is the largest multiple 
of 19 less than Ap +1. 


A B Cc 
10x 1— O=10 
10x 10— 95= 5 
10x 5— 38=12 
10x 12—114= 6 
10x 6— 57= 3 
10x 3— 19-π||] 
10x ll— 95=15 
10 x 15 — 133 = 17 

and so on. 
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The curiosity here is that the ‘units’ digits of column C are now 
those of column B in the last example whilst the ‘tens’ digits 
(which are to be ignored for the present purpose) occur in exactly 
the same positions as the remainders ‘1’ appeared in the previous 
column C. 

I am sure readers would find it entertaining to find for themselves 
the corresponding techniques for determining the recurring cycles 
of the fraction 1/29. There are, as one might expect, close relation- 
ships amongst all fractions of the form 1/(10n + 1) and it is not 
difficult, reasoning by analogy, to formulate general constructions. 


(2) These periodic sequences of digits relating to the common 
fractions have many curious properties, not all of them as predictable 
as one might expect and in order to observe some of them and to 
provide enough material for further experiment it will be useful to 


~have before us a sufficient number of examples. Table 3 lists the 


decimal equivalents of the reciprocals of the prime numbers from 
7 to 97 inclusive. 
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Table 3 
Recurrent Period of 1/p. (p—prime) 


7 142857 
1 -09 
13 -076923 
17 -0588235294117647 
19 .052631578947368421 
23 -0434782608695652173913 
29 -0344827586206896551724137931 
31 -032258064516129 
39) “27 
41 .02439 
43 -023255813953488372093 
47 -02127659574468085 10638297872340425531914893617 
53 0188679245283 
59 -016949152542372881355932203389830508474576271 1864406 
779661 
61 -016393442622950819672131147540983606557377049 1803278 
68852459 
67 -014925373134328358208955223880597 
71 = -01408450704225352112676056338028169 
73 = 01369863 
79 =-0126582278481 
83 -012048192771084337349397590361445783 13253 
89 = -01123595505617977528089887640449438202247191 
97 -0103092783505154639175257731958762886597938 144329896 
90721649434536082474226804 12371 1340206185567 


It will be seen at once that roughly half of these periods contain 
(p — 1) digits whilst the rest have some fraction of (p — 1). When 
the denominator of the original common factor is composite (say 
m) the number of digits in the period is, however, invariably less 
than m. For further reference and to provide material for any 
reader who wishes to pursue the subject, Tables 4 and 5 give the 
decimal equivalents of the reciprocals of some composite numbers. 
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361 


441 


Table 4 


Recurrent period of 1/m. (m—composite, ~ 2k, 5k, or k) 


‘025641 

‘019607843 1372549 
‘017543859649 122807 

‘015873 
‘0144927536231884057971 
‘012987 
‘0114942528735632183908045977 
‘010989 

‘010752688172043 


Table 5 
Recurrent period of 1/p?. (p 4 2k or 5k) 


020408 163265306122448979591836734693877551 
‘012345679 
‘0082644628099 173553719 
‘0059171597633136094674556213017751479289940828402366 
86390532544378698224852071 
‘0034602076124567474048442906574394463667820069204 I 52 
2491349480968858131487889273356401384083044982698969 
1937716262975785467 12802768 1660899653979238754325259 
5155709342560553633217993079584775086505 190311418685 
1211072664359861591695501730103806228373702422145328 
719723183391 
‘002770083 1024930747922437673130193905817174515235457 
0637119113573407202216066481994459833795013850415512 
46537396121 88365650969529085872576177285318559556786 
703601 1080332409972299 16897 5069252077562326869806094 
18282548476454293628808864265927977839335 18005540166 
20498614958448753462603878 1 16343490304709 14127423822 
71468144044321329639889 1966759 
‘002267573696145124716553287981859410430839 
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(3) With these tables we now have enough material to make a 
number of observations of the properties and methods of con- 
struction of these interesting periodic sequences. 

As we saw at the beginning of this chapter the fraction 1/19 is 
generated by multiplying successively by the “7 function (see 
Table 2) from the right and starting with unity. This method is 
quite generally applicable and is not restricted to prime denomina- 
tors. For instance the decimal equivalent of 1/39 as seen in Table 3 
is 0-025641, obtained by using the multiplier 4. 


Thus: 


4096 etc, 

.. . 1025641 
Multiplying this sequence by three gives . . . 3076923... (1/13). 
Of course the periodic part of the fraction 1/13 could have been 
obtained directly by multiplying as before by four but starting with 
the digit 3, thus: 


12288 etc. 
. . - 3076923 

Similarly, using x5, and starting with 7 we obtain the sequence 
. . . 142857 for the fraction 1/7. 

Using the above principles together with the ‘t’ functions in 
Table 2, we have then a general method for producing the periodic 
parts of any fraction whose denominator is not a multiple of two or 
five. 

There are many other ways of doing this. As an example of one 
employing a more orthodox arithmetic than the above we may 
note that on dividing seven into one hundred a remainder of two 
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is left; the following two hundred then has a remainder of four, 
and so on. The decimal equivalent of 1/7 can then be quickly 
written down by multiplying the first dividend by two, placing the 
result two places to the right and continuing in this way, thus: 

14 


448 and soon. 
1428571428... 
This example was chosen because it is often quoted as a curiosity ; 
in fact since 10 — 7 = 3 it is just as effective to use 3/10 as the 
multiplier. 


1428571... 

Being a straightforward substitute for long division there are no 
restrictions in this process and for example the fraction 1/8 can 
be calculated in an exactly similar manner. Here 8 = 10 — 2 and 
we have, 


. 1249999... -- .125 
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The method is obviously quite general, its main charm as a com- 
putational time-saver arising when the denominator of the required 
fraction or one of its multiples is close to and a little less than a power 
of ten. For instance, noting that 17 x 588 = 9996 = 10000 — 4, 
we can rapidly decimalise 1/17 by multiplying -0588 successively 
by 4/10000, thus: 


0588 
2352 
9408 
37632 
150528 
602112... 
05882352941176470588. . . 


The above examples lead back in a roundabout way to the second 
paragraph of this chapter, for clearly if a fraction’s denominator 
is a factor of 105 — 1, then the multiplier is unity. In such cases 
it is easy now to see that since, say, 999 = 37 x 27, then 1/37 = 
‘0270270. . . = 027, and as 99999 = 41 x 2439, then 1/41 = 
02439, 

The situation changes in a curious and less obvious manner when 
a denominator, or one of its multiples, is treated as exceeding a 
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power of ten. Taking again 1/17 as an example we have 17 x 6 = 
102, and we now use the multiplier 2/100 but this time alternately 
add and subtract the successive multiples. Thus: 


‘05882352 

ἣν 96 
‘0588235296 

- 192 
Ὁ58823529408 

+ 384 
‘05882352941184 

- 768 
0588235294117632 

Ῥ 1536 
‘058823529411764736 

- 3072 
‘058823529411764705... 


By a similar process, and noting that 7 x 143 = 1001, we have 
a still quicker method for determining 1/7, namely, 


4+ +143 

ἀρ: τ δὰ 
142857 

+ 143 


14285714... 


(4) It is in fact amazing how many algorithms can be devised for 
constructing the recurring periods of decimal fractions. Added to 
what has already been shown the following three examples will 
indicate something of the richness of this field. 
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Example 1. 1/31 


29 
2x 29+3 = 61 
2 x 61 + 29 = 151 
2 x 151+ 61 = 363 
2x 363+151 = 877 
2 Χ 877+ 363 = 2117 
2x 2117+877 = 5111 
2 x 5111 +2117 = 12339 
2 X 12339 + 5111 = 29789 
εν, 29032258064516129 
Example 2. 1/29 
31 
9 x 31 = 279 
9x279 = 2511 
9x2511 = 22599 
9 x 22599 = 203391 
9 x 203391 = 1830519 
. . . 551724137931 
Example 3. 1/67 
‘01 
49 
(1 + 49)/2 = 25 
(49 + 25)/2 = 37 
(25 + 37)/2 = 31 
(37 + 31)/2 = Se 
‘014925373134... 


(5) There are some still more remarkable ways in which these 
periodic sequences can be formed, illustrating in a striking manner 
how an unsuspected common link can sometimes be found between 
apparently unrelated branches of Number Theory. It will be more 
convenient, however, to examine these after introducing in the next 
chapter the curious additive series associated with the name of 
Fibonacci. 

From the constructions already given it will have become clear 
that fractions which have recurrent, or cyclic, periods can also be 
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expressed by the summation of infinite series. We have seen, in 
effect, that 1/7 (= 0-142857) is equivalent to 

A “fe ΕΣ + od + ha +... 

10 105. 10° 10: 
and also to 

Σ΄ Δ Se ἜΣ 

ΤΟΣ “+ To! + Τοῦ + ioe" 


which is the same as 


ἷ ss : ἢ + 

50° 502 (508 
and again, to 

13 143 143 


τοῦ τοῦ * 105 


It might also be noted that the series for 1/19 and 1/21 can be 
simplified into the following forms: 


1/19 = 1/20 + 1/202 + 1/20°+... 
1/21 = 1/20 — 1/202 + 1/208 -—.. . 


In general, if we express N as one of the forms (105 + a)/s 
the series for 1/N then becomes: 
1 5 sa sa sa 


(1) κεν - 


Ν τὸοΆι--α΄ 10. 1055 ἵ 1085 Τ᾿ 


(2) 1 4 Sa .C Ἐδη5 
Ν᾿ 10°+a 10" 10%" [085 "΄᾿ 
Some of the cyclic properties of recurrent decimals are of common 
knowledge, particularly the confined cycle which occurs when 1/N 
has a period consisting of N—1 digits. With a denominator of 
seven, for instance, the corresponding six fractions are: 


1/7 = 0-142857. . . 
2/7 = 0-285714. . . 
3/7 = 0-428571.. . 
4/7 = 0-571428. . . 
5/7 = 0714285... 


6/7 = 0:857142. . . 


| 
| 
| 
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This recycling of the same digits in the same order is to be found in 
the fractions 1/17, 1/19, 1/23, 1/29, 1/47, 1/59, 1/61, 1/67, 1/97,. . ., 
but it obviously cannot occur when the number of digits in the cycle 
is less than N — 1. Thus when N is 13 and 1/13 has the decimal 
equivalent of -076923 employing only (N — 1)/2 digits, we have 
two sets of recycling digits. 
1/13 = 076923. . ., 2/13 = -153846... 
3/13 = -230769. . ., 5/13 = -384615.. . 
4/13 = -307692.. ., 6/13 = -461538. . . 
9/13 = -692307. . ., 7/13 = :538461.. . 
10/13 = -769230. . ., 8/13 = -615384. . . 
12/13 = -923076. .., 11/13 = -846153.. . 


The two sets of numerators, (1, 3, 4,9, 10,12) and (2, 5, 6, 7, 8, 11), 
have interesting derivations and properties. For instance, if a list is 
made of the squares of the natural numbers 1, 2,3,. . ., and from 
these we subtract the largest multiple of thirteen that will leave a 
positive remainder, we have: 


1 4 9 16 25 36 49 64 81 
ee. SS DR RG 


re SS: 3 Bw Be 3 


which provides a recurrent cycle (in a different order) of the numbers 
1, 3, 4, 9, 10, 12, that is, the first set of numerators shown above. 

Again, the two sets, which we will call respectively A and B, 
have a curious internal property. Multiplying any two or more of 
the numbers of set A together and then subtracting a suitable 
multiple of thirteen leaves a remainder which is one of the numbers 
in the set. Thus (3 x 12) — 26 = 10, and (3 x 4 x 9) — (8 x 13) 
= 4, and so on. Moreover, on applying the same procedure to 
members of set B we again arrive at the numbers of set A, as 
follows (2 x 8) — 13 = 3. On the other hand if a number from 
set A is multiplied by one from set B the outcome is a member of 
set 8. 

It would be impracticable to pursue this line of thought further 
until some experience of congruences and quadratic residues has 
been gained but it must occasion some surprise that there is a direct 
link between recurrent decimals and the square numbers. 
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(6) And finally to round off this chapter I am going to revert to a 
construction method which the reader might care to experiment with 
and/or find an explanation. 

We see from Table 5 that 1/121 is equal to 


‘0082644628099173553719. . .. 


The following operations require no verbal description. 


8 == § 
8+18= 26 
26+18= 44 
44+ 18 τῷ 62 
62 + 18 = 80 
80 + 18 = 98 
98 +18 = 116 
116+ 18 = 134 
134+ 18 = 152 
152+ 18 = 170 
170 + 18 = 188 
188 + 18 = 206 
206 + 18 = 224 etc. 
826446280991735537190082.. . 


1+ 18= 19 
19 + 18 = 37 
37 + 18 = 55 
55 + 18 = 73 
73+ 18 = 9] 
91+ 18 = 109 

109 + 18 = 127 
127 + 18 = 145 
145 + 18 = 163 
163+18= 181 

181 + 18 = 199 etc. 


. « - 0082644628099173553719 


' 


— 
- 


OS eee naa 


-— eee - Gee-w eee 
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EXERCISES 
1, What are the decimal equivalents of 
3/41, 7/41, 13/41, 29/41, and 30/412 


2. Find empirically some algorithm which generates the recurrent 
period of the fraction 1/441. (= 1/21?), namely 


‘00226757369614512471655328798 1859410430839... 


This exercise might be expanded into a general search for 
algorithms among these periods. Anyone with a flair for empiri- 
cal observation will find ample scope for ingenuity in this 
field. 
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(1) One of the most exciting things about experimenting with num- 
bers is that one can so rapidly find oneself in strange and unfamiliar 
surroundings. And in such a situation there is little that can compete 
with the thrill of coming across some unexpected link with an 
apparently unrelated branch of mathematics. 

A particular and most interesting example of this is to be found in 
a study of the additive series originated by Leonardo de Pisa, 
nowadays better known as Fibonacci. This series, first defined in the 
early part of the thirteenth century, is usually presented in the form 
0, 1, 1, 2, 3, 5, 8, 13,. . ., each new term being the sum of the two 
preceding terms. In actual fact the fundamental properties of this 
important series would remain unchanged whatever integers were 
used to start it. Thus, observing the rule and giving any value we 
like to a, and ὦ, the series becomes: 

a, ὃ, (a + b), (a + 2b), (2a + 3b), (3a + 56), (Sa + 88), ete. 

For reasons which will be clarified later it will be found convenient 
to re-write the series in the form of two columns A and B, the 
respective terms of which are determined by the following instruc- 
tions: 

An +, = An + Bn, 
Bn +i, = Bn + An+1. 
Taking successive values from the above table and determining the 


decimal equivalents of the fractions Bn/An and (An+1)/Ba it will 
be found that these both converge towards the value 


1-6180339888 . . ., 
thus, 


89/55 = 161818. . ., 144/89 =1-61765... 
233/144 = 161806. .., 377/233 = 1-61798.. . 
610/377 = 1-61804. .., 987/610 = 1-61803. . . 


and so on. 


F 
| 
| 
) 
Ι 
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Table 6 
The Reconstructed Fibonacci Series 
n A B 
1 0 1 
2 1 2 
3 3 5 
4 8 13 
5 21 34 
6 55 89 
7 144 233 
8 377 610 
9 987 1597 
10 2584 4181 etc. 


We come now to a remarkable property of this series. The fractions 
An|Bn and Bn/(An+1) get closer and closer, as m increases, to 
061803... (for example 89/144 =0-61805...), while the 
fractions (An+1)/An and (Bn+1)/Bn approach 2-61803.. .. 

The Greeks knew this proportion (1°61803.. .) as the Golden 
Ratio and made frequent use of it in their architecture and sculpture. 
It seems to have a curiously aesthetic attraction and many famous 
canvasses have dimensions closely approximating to adjacent 
members of the series. The ratio is now more commonly referred 
to as the Golden Section and is variously indicated by the Greek 
symbols for ‘phi’ or ‘tau’, 

It is the ratio of a diagonal to a side in a regular pentagon and of 
the radius of a circle to the side of an inscribed regular decagon. 
It forms the ratio between the number of clockwise and anti- 
clockwise spirals in the seed patterns of such flowers as sunflowers 
and in those of the areoles of globular cacti. It can be used to generate 
the logarithmic spiral governing the growth of sea and snail shells 
and the distribution of leaves around the stems of many plants, and at 
least one dedicated research worker claims to have correlated it to 
the ratio of the height of women’s navels to their total stature. 

In spite, or perhaps because, of the richly coloured applications 
that have been claimed for this innocent looking series it has come 
to be regarded more as a mathematical curiosity than as a subject 
for serious study and has received scant attention in most writings 
on the higher arithmetic. In what follows I hope to modify this 
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impression and to show that there is a sound mathematical content 
in series of this form. 


(2) First to deal with the fairly well-known properties of the Golden 
Ratio. Partly to avoid conflict with the adherents of ‘phi’ and ‘tau’, 
and partly to cope with the limitations of my typewriter, I shall use 
‘F’ to indicate 1-61803. . .. 

It is readily verified that (1-618. . .)5 =2°618..., and that 
1/1618... = 0-618. . .. Thus F is a number whose reciprocal 
is obtained by subtracting unity and which is squared by adding 
unity. Both of these statements, 1/F =(F—1) and FP? =F+1 
produce the quadratic equation 

F?— F—1=0 
Solving this by normal algebraic procedure we obtain: 


= en 
ae) = (1 + +/5)/2 = 1-618. . .. 
It follows that 1/F = (4/5 — 1)/2 = 0-618. . ., and since +/5 is 
irrational, that is its decimal part neither terminates nor recurs, then 
F is also irrational. 

Both F and the series from which it is derived have strange pro- 
perties. For instance, keeping one eye on Table 6, note the following 
relationships: 


F= 


2x 5— #2= 2?— 1x 3=1 
5x13— 883. 5S?— 3x 8=1 
13 x 34 -- 213 τῷ 133 -- 8x21=1 
34 x 89 — 552 = 342 — 21 x 55-Ξ1 and so on. 


The reader may well find it interesting to search for other associations, 
perhaps using different starting numbers for ‘A’ and ‘B’, but I 
think it is safe to predict that he will find it difficult wholly to escape 
from the sequential integers of the series. 

Even the ascending powers of F cannot disengage themselves as 
will be seen from the following relationships, 


F= F 

F= F+1 
Fe= 2F+1 
Fi= 3F+2 
F§= 5F+3 
Fe= 8F+5 


Fi =13F +8 and so on. 
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But perhaps the most unexpected property of this series emerges 
when the two numbers 1-61803... and 261803. . . are multiplied 
in turn by the natural numbers 1, 2, 3, 4,. . ., the fractional parts 
of the products being ignored. The first few pairs of numbers formed 
in this way are: 


(1, 2), (3, 5), (4, 7), (6, 10), (8, 13), (9, 15), (11, 18), 


and it can be shown that however far the operation is carried no 
integer is repeated and none omitted. 


(3) References to the ‘magical’ properties of the Golden Section are 
often to be found in articles dealing with mathematical puzzles and 
paradoxes and it is often claimed, or at least implied, that it occupies 
an unique position in regard to the peculiarities which have just 
been mentioned. This may be true in some senses but let us return 
to our presentation of the Fibonacci series in two columns and try a 
slight change in its formulation. 

The ‘A’ and ‘B’ columns will now be constructed according to 
the rules: 


Aa+1 = An + Bn 
Bn+1 = An + An+1 


We then have, 
Table 7 

n A B 

1 0 1 
2 1 1 
3 2 3 
4 5 7 
5 12 17 
6 29 41 
7 70 99 
8 169 239 
9 408 577 
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One of the advantages of using a two column presentation of such 
series will now become apparent; stretched out in a single line as 
Fibonacci’s series is usually written it would be difficult to recognise 
this one as purely additive. However by its method of construction 
it most certainly is, and from an arithmetical point of view it is 
equally interesting. 

Applying to the values in Table 7, the technique which provided 
the F number 1-61803... and its derivatives we find that the 
ratio Bn/An approaches, as ἡ increases, the limit 141421356237 

. «, Which is of course the decimal expression of 4/2. 

Furthermore the ratios (An+1)/An and (Bn+1)/Bn, as n in- 
creases, approach the value +/2 + 1, (Bni1)/An gets nearer and 
nearer to 1/2 +2, and so on. Incidentally, these ratios and the 
general properties of this series are, as in the Fibonacci construc- 
tion, quite unaffected by the choice of the two starting numbers. 

This latter series also shares what seemed to be a property peculiar 
to the ‘F’ numbers. Multiplying the values 1:4142... and 
3-4142. . . successively by the whole numbers 1, 2, 3,4,... and 
discarding the fractional parts of the products we again find a 
sequence of pairs of numbers which includes all the natural numbers 
without duplication. The first few of these pairs are: 


(1, 3), (2, 6), (4, 10), (5, 13), (7, 17), (8, 20), (9, 23),. . . 


and here we see that the differences between the numbers in brackets 
are successively 2, 4, 6, 8, etc. 


(4) Striking as these properties and the relationships between the 
two series are, there is a still more remarkable connection to be 
noticed between them and some of the periods of the recurrent 
decimals examined in Chap. 3. Four examples from my own 
observations will be given and I have no doubt the reader will find 
this an exciting field for further research. 

(A) If we take the original Fibonacci series, 0, 1, 1, 2, 3, 5, etc., 
and divide each term in turn by 10, 100, 1000, and so on, and then 
sum the results we have: 
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‘0112358 


‘011235955056179 . . . = 1/89 (see Table 3). 


(B) If, on the other hand, the alternate terms of this series—as shown 
in column A of Table 6—are treated in the above manner we 
arrive at the following summation: 


‘0138 


14930352 
39088169 
102334155 


014084507042... . = 1/71. 


(C) Now we stated earlier—and this is easily demonstrated—that 
the ratio between successive terms of the Fibonacci series approaches 
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F whatever integers are used to start the sequence. Again using the 
above process but this time starting with the numbers 5, and 26, 
and dividing the successive terms by 10?, 104, 10° etc., we now have: 


0526315788 
145 


‘05263157894736842105263. . . = 1/19 


(D) Turning now to (4, 3), in which a series producing the ratio 
141421. . . (4/2) was described, it will be seen that this sequence 
of numbers also has its surprises. Taking alternate terms—that is the 
consecutive lines of column A in Table 7—and treating them as in 
(A) and (B) above we obtain the decimal equivalent of the fraction 
1/79, thus: 


13860 and continuing, 
33461 
80782 
195025 
470832 
1136689 
2744210 


‘01265822784810. . . = 1/79. 


, 


- --- 


- gt eee 
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I think this link with the recurrent periods of ‘prime denominator’ 
fractions presents additive series in a new light and certainly there is 
scope here for some research. 


EXERCISES 
1. Calculate 6F?/5 to five significant figures. 
(F = 1-61803. . .) 
2. What is the sum of the first m Fibonacci numbers? 


5 ONE ONE ONE ....... ONE 


(1) In the previous chapters we have frequently run up against those 
numbers which belong to one or the other of the forms 10° — 1 
or 10+ 1. Expressed in digital language these are respectively 
of the general types 


9999. . .99 and 1000. . . Ol. 


We have seen that rapid checks of divisibility can easily be applied 
to the factors of such numbers and that they have an important 
bearing upon the structure of recurring decimals; for these and a 
variety of other reasons it will be of interest to examine them in 
greater detail and to discover some of their special properties and 
behaviour. 

All integers of the form 415 —1 are divisible by A —1 and 
therefore in dealing with numbers of the form 10" — 1 it will be 
most convenient first to remove the factor 9 and concentrate upon 
the sequence 1, 11, 111,..., 111. ..1. In what follows such 
integers will be referred to as ‘J’ numbers, and in particular J, 
will be used to indicate (105 — 1)/9, or the integer 111...1 
having n digits. Thus, for example, J5 = 11111 = (10° — 1)/9. 

It will also be convenient to use a similar convention to describe 
integers of the form 10°+ 1; these will be referred to as ‘J’ 
numbers, the subscript again being equal to n. Thus Js, or 10° + 1 
represents 100001, in this case being the number of zeros plus one. 

By definition J, Χ Jn = Jan and in general it is easily seen that 
when n is composite J, can always be resolved into at least two 
factors, as for example: 


112 = 11 x 10101010101 = 111 x 1001001001 
= 1111 x 100010001 = 111111 x 1000001 


In can therefore only be prime when 7 is prime and indeed it 
is one of the most remarkable properties of both the ‘J’ and ‘J’ 
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numbers that they seem to be extraordinarily free from primes. 
(I say ‘seem’ because the existing tests of primality become exceed- 
ingly difficult to apply when » exceeds 37.) 

Iz (= 11) is certainly prime, Ji9 is generally accepted to be prime, 
although I have not found any reference to a positive proof of this. 
On the other hand J23 is still widely quoted as being prime in 
spite of a somewhat dubious ancestry and a proof of its composite 
character by D. H. Lehmer in 1927. Beyond this there are no more 
prime J,’s up to, at least n = 137. 


(2) One might expect that numbers constructed in such an orderly 
and regular fashion would exhibit some sort of pattern in the dis- 
tribution of their factors. As we shall see later this is to a certain 
extent true but at first sight complete factorisation presents a 
chaotic picture. 

The following Tables 8 and 9, list the factors of the smaller 7 and 
J numbers. I think these factors are all prime but am not completely 
certain about the large factors of Jes, and those of Jy, where 
n = 20,22, and 23. Perhaps after studying the factorisation methods 
which will be described later some readers may find the subject 
interesting enough to try their hand at cracking these somewhat 
formidable numbers. 
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Table 8 


Factors of J, (- εὐκάμέξη =) 


Factors 
11 (prime) 
3.37 
11.101 
41.271 
3.9044 13.37 
239.4649 
11.73.101.137 
3.3.37.333667 
11.41.271.9091 
21649 . 513239 
3,7.11,13.37.101.9901 
53.79 .265371653 
11.239.4649 .909091 
3.31.37.41.271.2906161 
11.17.73.101.137.5882353 
2071723 . 5363222357 
3.3.7.11.13.19.37.52579 .333667 
1111111111111111111 (thought to be prime) 
11.41.101.271.9091 .99009901 
3.37.43 .239. 1933 .4649 . 10838689 
11.11.23.4093.8779.21649 .513239 
(Composite but factors not known) 
3.7.11.13.37.73.101.137.9901 .99990001 
41 .271 .100001000010000100001 
11.53.79.859 .265371653 . 1058313049 
3.3.3.37.757.333667 .440334654777631 
11.29.101.239.281.4649.909091 . 121499449 
3191 . 16763 .43037 .62003 . 77843839397 
3.7.11.13.31.37.41.211.241.271.2161 .9091 .2906161 


11.17.73.101.137.353.449 641.1409 .69857 .5882353 
11.103 .4013 .2071723 .5363222357 .21993833369 
3.3.7.11.13.19.37.101 .9901 .52579 .333667 .999999000001 
11.1111111111111111111.909090909090909091 


a “«-........-... “Ὁ 
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Incidentally, J3g is unique among the numbers listed above in that 
it is the only one, with even, that has three prime factors. 


Table 9 
Factors of J, (= 10" + 1) 
Factors 


11 (prime) 

101 (prime) 

7.11.13 

73.137 

11.9091 

101.9901 

11.909091 

17.5882353 

7.11.13.19.52579 

10 101.99009901 

11 = 11.23.4093.8779 

12 73.137.99990001 

13 11.859. 1058313049 

14 29.101 .281.121499449 

15 -7.11.13.211.241.2161.9091 

16 353.449.641.1409.69857 

17 11.103.4013.21993833369 

18 101.9901 .999999000001 

19 11.909090909090909091 

20 73.137.9999000099990001 

21 = 7.7.11.13.127.2689.459691 .909091 
22 89.101.1112470797641561909 

23 -11.47.139.2531.549797184491917 
24  17.5882353.9999999900000001 

25 :11.251.5051.9091 .78875943472201 


(3) Of course, if we are content to express these numbers by, at 
most, only two factors a much more regular pattern emerges, as 
exemplified in the following breakdown of the J numbers. 


won aA Un hh W ND μὰ 


J2 101 J3 11.91 

Ja 10001 Js 11.9091 

οἷς 101.9901 .Ἶ) 11.909091 
Js 100000001 Jo 11.90909091 
Jio 101.99009901 etc. 
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In fact these numbers containing only the digits 0, 1, and 9, 
arranged in various patterns offer in themselves an interesting sub- 
ject for research; one set in particular provides an excellent example 
of the dangers of drawing hasty conclusions in enquiries of this kind. 
Taking in sequence one of the factors of Jgn (or Jen) we have: 


J3 91 composite 
Je 9901 prime 
Jo 999001 composite 
etc. 99990001 prime 


9999900001 composite 
999999000001 prime 

99999990000001 composite 
9999999900000001 prime 

999999999000000001 composite 


One might be tempted to infer that the next number in the sequence 
(99999999990000000001) is prime but fortunately as it must be a 
factor of Jo it is not difficult to prove that it is composite. 

Since 60 contains the factors 2, 3, 4, 5, 6, 10, 12, 15, 20, 30, it 
follows that Jgq must contain all the factors appearing in the 
numbers Je, Js, J4,. . . 180, namely, (11), (3.37), (101), (41.271), 
(7), (9091), (9901), (31.2906161), (99009901), and (211.241.2161). 
It is easy to see that if these factors are multiplied together the 
product must have more than forty digits and hence Jg9 cannot have 
a single prime factor of twenty digits. 

The number 99999999990000000001 must therefore have at least 
one of the above primes as a factor. 


(4) It will be seen from Tables 8 and 9, or for that matter readily 
deduced from first principles, that divisors of J, numbers are also 
divisors of Jmn where m is one of the natural numbers 1, 2,3,.. .. 
Consequently when dealing with the J and J numbers we need 
really only be concerned with those primes which appear as divisors 
for the first time. These will recur in all multiples of m and further 
research on numbers of this form will be made easier by abstracting 
the primary factors and presenting them in tabular form as in the 
following Table 10. 
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Table 10 
Primary Factors of ἴῃ 


Factors 


21649 . 513239 

9901 

53.79 .265371653 
909091 

31.2906161 

17. 5882353 
2071723 . 5363222357 
1952579 
1111111111111111111 
99009901 

43 .1933 . 10838689 

23 .4093 .8779 

not known 

99990001 
100001000010000100001 
859. 1058313049 

757 .440334654777631 
29.281. 121499449 
3191. 16763 .43037 .62003 . 77843839397 
211.241.2161 

2791.N 
353.449 .641 . 1409 .69857 
67.N 

103 .4013 .21993833369 
71.N 

999999000001 
2028119.N 
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This table provides a number of working examples of one of the 
most important theorems in the Theory of Numbers and at the same 
time an insight into a puzzling feature of the varying periods of 
recurring decimals. 

In the first place, remembering that factors of J, are also factors 
of Inm we see that: 


7 is a factor of Jg 
11 is a factor of [πὸ 
13 is a factor of Jie 
17 is a factor of Ji¢ 
and so on. 


And generally, when 7 is prime it is a divisor of Jn-1. 

Bearing in mind how the J, numbers are constructed this is seen 
to be a particular example of ‘Fermat’s Theorem’ which states that 
if p is prime and W is not a multiple of p, then N?-+—1 is 
divisible by p exactly. This theorem is of fundamental utility in the 
exploration of numbers; proofs, depending usually either on binomial 
expansion or congruence theory, will not be given here as they are 
to be found in almost every textbook on Algebra. 

The table also serves to illustrate the following lesser—but 
still important—theorems. Here again the proofs, which though not 
difficult depend upon a sequence of lemmas, are left to the formal 
textbooks. (For instance, Advanced Algebra, Barnard and Child 
Macmillan & Co. Ltd.) 


(1) If n is an odd prime, every prime factor of N” — 1 which is not 
a divisor of N — 1 is of the form 2kn + 1. 


(2) If p is a prime factor of N” + 1 it is also a factor of Νά - 1 
where d is the greatest common divisor of m and }(p — 1). 


(3) If nm is a prime, every odd prime factor of N“ + 1 which is not 
a divisor of N + 1 is of the form 2kn + 1. 

(4) If πὶ is a prime, every prime factor of N®” — 1 which is not a 
divisor of N** — 1 is of the form kn? + 1. 


(5) Every prime factor of N2” + 1 is of the form 2"+1. k + 1. 
ry p 


These theorems are of value in limiting the choice of possible fac- 
tors of both J and J numbers and as we shall see later it is possible 
to narrow the field still further. 
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(5) Before pursuing this line however, it will be of interest to note a 
connection between recurring decimals and Table 10. Referring 
back to Table 3 and its following paragraph we recall that the 
recurrent period of the decimal form of 1/p contains a number 
of digits which is equal either to p—1 or to some divisor 
of p—1. Thus 1/7, (= -142857) has six digits whilst 1/53, 
(= 0188679245283) has 52/4 = 13 digits. At first sight the period 
lengths appear to be decided in an entirely random manner. 

But the factors of J, are in effect the factors of 999. ..9 
(having πὶ 9’s) and the periods of recurring decimals are exact 
divisors of numbers of this kind. In consequence the reciprocals of 
the factors shown in Table 10 have the corresponding m number of 
digits in their recurring period. Thus both 1/41 (= -02439), and 
1/271 (= -00369) have five-digit cycles, 1/7, and 1/13, six-digits, 
and so on. Furthermore in the whole infinity of the natural numbers 
there can only be one whose reciprocal has exactly 2, 4, 10,... 
19,. . . 36, etc., digits. 


(6) The factorisation of J numbers presents a challenge which is 
difficult to resist and although we shall deal with general methods of 
factorisation later this is an appropriate point to mention some aids 
which can be applied to these numbers in particular. Unfortunately 
it is outside the scope of this book to discuss the theory of Quadratic 
Residues (that is, the remainders left after subtracting a given prime 
or its multiples from the sequence of square numbers) but one of its 
consequences must be mentioned in this context. 

Numbers which have a Quadratic Residue of 10 are of the general 
form 40k + 1, 3, 9, 13, 27, 31, 37, 39, whilst those with a Q.R. of 
—10 are of the form 40k + 7, 11, 17, 19, 21, 23, 29, 33. It can be 
shown that primary factors of J, (1.6. m odd) must be of the first 
form whilst those of J, are of the second. This, I am afraid, will 
have to be taken on trust as far as formal proof is concerned but the 
factors given in Table 10 and a subsequent table (12) will show a 
practical confirmation of these statements. 

The important thing is that we now have two independent ‘forms’ 
of factors of these numbers and these can be combined to give a 
much more economical approach to their factorisation. 

To take a specific example let us examine the possible factors of 
J; (= 1,111,111). These must be of one or more of the forms 40s + 
1, 3, 9, 13, etc., and at the same time of the form 1417 + 1: equating 
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these forms in turn we have: 


(A) 14¢+1=40s+1, and hence 7t = 20s. This equation is 
satisfied when ¢ = 20,40,. . . 20k. Substituting, wehave 14(20k) + 
1 = 280k + 1. 


(B) 14¢+1=40s+3, thus 71 Ξε 205 1, which is satisfied 
when ¢ = 3, 23, 43,. . .20k - 3. Then 14(20k + 3) + 1 = 280k 
+ 43. 


Continuing in this way it is not difficult to determine that the 
factors of J7 must be of one of the forms 


280k + 1, 43, 71, 169, 197, 239, 253, 267. 


There are 168 odd primes less than the square root of J7 but 
by using this elimination process one factor (239) is disclosed 
at the fourth trial. The other factor (4649) is equal to (280 x 16) + 
169. 

The following table has been constructed on these lines and should 
prove useful to anyone intending to pursue the subject. 
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Table 11 


‘Forms’ of Factors of J, 


n 
5 40k+1, 31. 
7 2ae+1, 4%, Ti, 168, i197, 239, 253, 267. 
> SORE, 57, 1688). 199; “253; 2h, 269 1 
ll 440%+1, 67, 85. I, 133, 199, 243, 397. 
Ὁ fer. 2 3 Ὁ Tf 2a as Bet. 
155. 120kK+1, 31. 
17 680k +1, 239, 307, 409, 443, 477, Sil, 613. 
Sy etl FF TT, TS, Bt 55, Ge Tex 
21 6M +1, 45, 100, 253, 7, G31, 755, 5, 
23 Fieri, ΤΣ 231,. 277, 33, ea, Bi, Bes. 
25 200K+1, 151. 
27 1080kK+1, 163, 271, 649, 757, 919, 973, 1027. 
29 1160K +1, 117, 523, 639, 813, 871, 929, 987. 
31 1240kK+1, 187, 249, 311, 373, 559, 683, 1117. 
oS Ise +1, δὴ, 133, 199, 37, Ὁ. SA, 1H 
35 280K+1, 71. 
37 1480k+1, 519, 667, 889, 963, 1037, 1111, 1333, 
39 1560K+1, 79, 157, 391, 547, 1093, 1249, 1483. 
41 1640k+1, 83, 329, 493, 1067, 1231, 1477, 1559. 
43 1720Κ -ἰ 1, 173, 431, 517, 603, 689, 947, 1119. 
45 360k+1, 271. 
47 1880k +1, 283, 471, 1129, 1317, 1599, 1693, 1787. 
49 1960k +1, 197, 883, 1079, 1373, 1471, 1569, 1667. 


(7) Even with the aids described above, however, the factorisation of 
I, and Jy becomes extremely difficult as m increases and there is 
scope here for further research into links between possible prime 
factors and n. As a start in this direction it might be worth while 
to note that although J, is not always composite it is obvious that 
every prime is a factor of some I, (see the opening remarks on 
recurring decimals in Chap. 3). It may therefore be of interest to 
approach the problem, empirically, from the rear. 

In the following table the odd primes up to 1051 are listed, 
together with the » values of J, for which they are factors. Thus, 
for example, it will be seen that when p = 41, the corresponding 


Table 12 
Values of πα for which J, is a multiple of p 
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value of n is 5, indicating that 41 is an exact divisor of Js, or 
11,111. 

(It is worth noting that when p is of one of the forms 40k + 7, 
11, 17, 19, 21, 23, 29, 33, then m is always an even number and 
consequently p is also a factor of Jn.) 


(A still more important observation can be divided into two parts: 


(a) When p is of the above form, 40k + 7, 11, 17, etc.—marked 
thus * in the table—then n is equal to (p — 1) divided by an odd 
number. For example, (17 — 1)/1 = 16; (73 — 1)/9 = 8, etc. 

(δὴ) When p is of the form 40k + 1, 3, 9, 13, 27, 31, 37, 39, n is 
equal to (p — 1) divided by an even number, as (31 — 1)/2 = 15; 
(173 — 1)/4 = 43.) 


Note also: p = 493121, n= 67 (Brillhart) 
497867 89 (Brillhart) 

18797 127 (McCullough) 

80173 131 (McCullough) 

2467 137 (McCullough) 

12517 149 (McCullough) 

12671 181 (McCullough) 

52009 197 (McCullough) 
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(1) We saw in the last chapter that when ἡ is composite it is always 
a simple matter to find at least two factors of J, and although it may 
not be easy to say whether these are prime or composite the process 
does reduce the magnitude of the numbers to be tested. The examples 
given employed only the digits 0,1 and 9, but other factor patterns 
can be found when 7 is of certain specific ‘forms’. 

For instance, we have: 


(A) When n = 3k, 


Ig = 3.37 
Ie = 33.3367 
Ig = 333.333667 


Ihe = 3333.33336667 and so on. 


(B) When n = 6k, 


Ie = 91.1221 

Tig = 9901.11222211 

Tyg = 999001.111222222111 

Ieq = 99990001.1111222222221111 ete. 


(C) When n = 4k + 2, 


Ig = Ill 

Ie = 91.1221 

Iho = 9091.122221 

14 = 909091.12222221 

Tig = 90909091.1222222221 etc. 
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Associated with the above the following curiously regular ‘irregu- 

larities’ are worth noting: 

3367 = 37.91 

33336667 = 37.900991 

3333366667 = 37.90090991 

33333336666667 = 37.900900990991 

3333333366666667 = 37.9009009099099 1 

33333333336666666667 = 37.90090090099099099 1 

3333333333366666666667 = 37.90090090090990990991 etc. 


A more intriguing partition of factors of certain J, numbers is to 
be found in a system which is not completely consistent. When 
n = 5k, In can be divided into the following pairs of factors: 


5 41.271 
10 451.2463661 
15 4551.24414658561 
20 45551.243926831707561 
25 455551.2439048780731709756 + 55555 


30 4555551.24390268292707317097561 

35 45555551.243902463414636585366097561 

40 =455555551.2439024414634146585365856097561 

45 = 4555555551.24390243926829268317073170756097561 

50 = 45555555551.2439024390487 804878073 17073170975609756 + 
5555555555 


and so on, with every J25x failing to conform to the pattern. 
(Note that 455551 = 412. 271.) 

This peculiarity is not confined to members of the 755 clan as 
can be seen by making a similar partition when a = 7. 


7 239.4649 
14 2629.4226364059 
21 26529.4188288707117159 
28 265529.4184518870297071548159 
35 2655529.4184142259832640167405858159 
42 26555529.4184104602514644355648539748958159 
49 265555529.4184100836820125523016736402092050251046 +- 
7777777 


and so on. Again, 265555529 = 2392 . 4649. 
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(2) Similar investigations into the construction of the J numbers 
must, I am sure, be equally rewarding and I recommend this as a 
project for the reader. In case however, the large numbers with 
which we have been dealing have been a little over-powering this 
chapter will be rounded off with a few observations which do not 
demand quite so much muscular mathematics. . 

It may seem hardly worth remarking that since 41, for instance, 
is a divisor of 11111, it must also divide 4111111, but it is perhaps 
not quite so obvious that it is also a factor of 4000001, 400000000001, 
etc. Similarly 37 is an exact divisor of 31117, 31111117,. . . and 
30007, 30000007, etc. 

The reasons for these ‘phenomena’ are easily discovered but what 
follows is not so obvious. Since, for example, 4649 is a factor of 
Ty it is also an exact divisor of 40000000649, 46000000049, 
46400000009 and all other members of the same pattern containing 
7k zeros (or for that matter 7k ‘ones’). 


(3) We have seen that the J and J numbers can present problems 
of interest and considerable complexity and that the determination 
of their prime or composite character is in general extremely diffi- 
cult; on the other hand their construction by additive methods is, 
as might be expected, not only simple but open to a variety of ap- 
proaches. Some of these relating to the 7 numbers are given below; 
the reader may find it diverting to extend these examples or perhaps 
to find similar expressions for the J numbers. 


(4) 1=(6—S/l=(7—4/3 =( — 3/5 =O —2)/7... 
11 = (62 — 52}}1 = (72 — 42/3 =... 
111 = (562 — 55%)/1 = (572 — 54/3 =... 
1111 = (5562 — 555%)/1 = (5572 — 5542)/3 =... and soon. 


(B) 11 = (105 — 1)/9 
111 = (602 — 512)/9 
1111 = (5602 — 5512)/9 
11111 = (55602 — 55512)/9 
111111 = (555602 — 555512)/9 etc. 
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(C) 11=(9 x 9) + 7)/8 
111 = ((98 x 9) + 6)/8 
1111 = ((987 x 9) + 5)/8 
11111 = ((9876 x 9) + 4)/8 


111111111 = ((98765432 x 9) + 0)/8 
1111111111 = ((987654321 x 9) — 1)/8 


(D) 11=1x9+2 
111 =12x9+3 
1111 = 123 x 9+4 
11111 = 1234 x 9+ 5 


111111111 = 12345678 x 9+ 9 
1111111111 = 123456789 x 9 + 10 


(FE) 11- 63 -- 58 τι 42-5 = 3242 
1111 = 562 — 452 = 342 — 45 = 332 + 22 
111111 = 5562 — 4452 = 3342 — 445 = 3332 + 222 
11111111 = 55562 — 44452 = 33342 — 4445 = 33332 + 2222 etc. 


(F) 11=62+75=72+ 62 
11111 = 662 + 6755 = 672 + 6622 
1111111 = 6662 + 667555 = 6672 + 666222 
111111111 = 66662 + 66675555 = 66672 + 66662222 etc. 


(G wW=(3x4-1 
1111 = (33 x 34) —11 
111111 = (333 x 334) — 111 
11111111 = (3333 x 3334) -- 1111 ete. 


(H) 1=1 
11=2+49 
111=3+(3 x 9)+92 
1111 =4+4+(6 x 9)+(4 x 92) - 98 
11111 = 5+ (10 x 9) + (10 x 92) + (5 x 98) +94 ete. 


Here, in effect, J, has been expressed in the scale of nine, and it 
will be seen that the coefficients are those of the binomial expansion 
of (1 + x)". 

If should be noted—though I hope without surprise—that if the 
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sequences of numbers defined by 2 + 1 are expressed in the scale 
of 2, they become: 


Qn ἢ Qn 4] 
n S10 Se Sio Se 
1 1 1 3 11 
2 3 11 5 101 
3 7 111 9 1001 
4 15 1111 17 10001 etc. 


And finally, although it is outside the intended scope of this book 
I cannot resist adding this little gem for the entertainment of those 
readers who are familiar with the operation of determinants: 
In In+1 In+2 = 0 
In+1 Inte In+ 
In+2 In+3  In+ 


In expanded, and far less picturesque form, this becomes (for 
compactness putting 0, p,q, 7, in the place of n + 1, + 2, etc.): 
InIp X Tp — Iq?) + To(Ig X Ip— Tr X 10) + Ip(lo X Ig — Ip”) = 0 

Testing this for n = 1, for example, we then have: 

(111.11111 — 11112) + 11(1111.111 — 11111.11) + 
111(11.1111 — 111%) = (1233321 — 1234321) + 
11(123321 — 122221) + 111(12221 — 12321) = 


—1000 + 11(1100) + 111(—100) = —1000 + 
12100 — 11100 = 0. 


EXERCISES 


1. What factors are common to both 1001001001 and 100010001? 
2. Why must 90909091 be composite? 
3. What are the factors of 33336667? 


7 PRIME NUMBERS 


(1) Primes, those integers which have no factors other than them- 
selves and unity, have fascinated and tortured amateurs and some 
of the greatest mathematical minds in history alike for centuries. 
And the total of our knowledge of them is still disappointingly 
small. 

We know that the primes continue indefinitely—that there is no 
‘largest’ prime; we know that no algebraic expression can be written 
that generates only primes. We know an analytic formula which will 
tell us how many primes to expect less than any given number, but 
only to a very rough approximation, and a procedure giving an 
exact figure at the expense of extremely laborious computations. 
We know of some expressions in x (variable) that contain in- 
finitely many primes but these are mostly trivial and though we 
suspect others they have not yet been proved. 

A few tests which distinguish between prime and composite 
numbers are known but they are difficult to apply to integers of even 
moderate size. And, broadly speaking, that’s about the lot. 

Oddly enough, in spite of what has been said, we can compose 
strings of consecutive integers which are composite for as many terms 
as we like to specify although to be frank this is of little value in the 
theory of primes. 

Having said this let us look at some of the intriguing ways in 
which the problem of the primes has been attacked. 

Proof of the infinity of the primes dates back at least to Euclid 
whose reasoning is a model for all time of succinct elegance. Suppose, 
he argued, that P is the largest prime: then if P and all the smaller 
primes are multiplied together and 1 is added to the product, then 
this new number cannot be exactly divisible by P or any lesser 
prime. It must therefore either be itself a prime or divisible by some 
prime greater than P; in either case there is a prime greater than P. 

Although there can be no formula which produces only primes— 


Prime Numbers 69 


proof of this may be found in any college algebra—there are some 
expressions which have an almost magical facility up to a point. 
For instance, the sequence p + 0, p + 2,p + 6,. . . (which is the 
same thing as p + n? — n,n taking the successive values 1,2, 3,.. 3) 
obviously becomes equal to n® when n =p but up to this point, 
when p = 3, 5, 11,17, and 41 we have: 


$7 8 TFT Φ Mi. tits 
511 7 Se 955. ..Σ 
7. 5. 2 & 421 Ἔσο 
> 31 2 GG 41 1572 
41 47 71 503 1447 
53 59 83 547 1523 .-.393... 2 
67 73 97 593 1601 = 40? +1 
83 89 113 641 1681 = 41? 


149 173 797 
is i 809 
199 223 911 
ant 253 $F 


257 281 1033 
289 313 1097 


All these integers except the terminal ones are prime. It is diffi- 
cult to believe that it is only by chance that these sequences pick 
their steps so daintily without once treading on a composite number 
but so far no one has been able to find a prime greater than 41 which 
produces a further series of this kind. 

Other expressions are occasionally ‘discovered’, but these arise 
from the fact that nm can be given negative values or be replaced by, 
for example (mn — 40), and in this way it is possible to devise a 
confusing array of ‘different’ formulae. Thus, for instance, we can 
write: 


n2 + 3n + 43 
n? + 5n + 47 
n2 + 7n +53 and so on, 
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or "3 — 79n + 1601 
2 
or n? +- an + a) (when a is odd). 


All these formulae produce exactly the same sequence of prime 
numbers as the original. 


(2) A fundamental property of prime numbers is conveyed in the 
statement, ‘If m is any integer and p is a prime, then n? —n 
is exactly divisible by p.’ 

Over 2000 years ago the Chinese knew this to be true in the special 
case of n= 2, but nowadays the above is known as ‘Fermat's 
(lesser) Theorem’. Fermat first declared the general form—without 
proof—in a letter to a fellow mathematician in October 1640. 
Removing the factor n we obtain the form in which the theorem is 
more commonly stated, namely; ‘If p is a prime number, and n 
an integer not a multiple of p, then n?-! — 1 is exactly divisible 
by p.’ 

The first formal proof of this theorem was given by Leibniz 
(1646-1716). In an age in which algebraic symbolism was in a 
primitive state and which consequently required propositions of 
this sort to be argued out verbally the proof must have presented 
formidable difficulty: today it is disposed of in a few lines in any 
school algebra. 

The theorem permiates the Theory of Numbers like a mycelium 
but unfortunately its converse cannot be relied upon as a test of 
primality without some highly sophisticated refinements. 

Another important theorem connected with prime numbers is that 
known as ‘Wilson’s’. This, reminiscent of Euclid’s proof of the 
infinity of primes, states that if, and only if, p is a prime, then 


C945... p= D+ 


is exactly divisible by p. (2.3.4.5... .m is known as ‘factorial π᾿ 
and is written n!). In most Algebra schoolbooks the theorem is 
stated thus: 

If p is a prime number, then (p — 1)! +1 is divisible by p. 
Contrary to Fermat’s this theorem is only true when p is prime 
as can easily be shown. For if p has a factor, say m, then m is 
less than p and must therefore divide (p — 1)!. The same goes for 
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its other factor and hence neither of them can be divisors of (p — 1)! 
+ 1. 

Factorial numbers increase very rapidly even for quite small 
values of m, as can be seen from the table in the Appendix, and 
consequently Wilson’s theorem has no value at present as a test of 
primality. On the other hand it has some interesting consequences. 
For instance it is not difficult to derive from it that if p and p +2 
are both prime then 


2(p — 1)! + 1 is divisible by p + 2. 


Anything bearing on ‘pairs’ of primes such as 5/7, 11/13, 17/19, 
59/61, 115301/115303, 100004561/100004563,  1000000009649/ 
1000000009651, etc., is worth noting since the theory of these 
numbers has scarcely been scratched and it is not even known 
whether they continue indefinitely. 

Perhaps if some technique for dealing with large factorials can be 
developed, comparable say with that now available for use against 
the converse of Fermat’s theorem, Wilson’s might yet disclose new 
horizons. 

Although this property of factorials does not appear to provide an 
effective weapon against the primes factorials do enable us to make 
the most fantastic demands upon composite numbers. Lists of 
prime numbers display irregularly spaced gaps indicating strings of 
consecutive composite integers; these appear to be randomly dis- 
tributed and increase in length as the numbers become larger. 
Nevertheless it would require a long list of primes and quite a lab- 
orious search before we could find say, fifty numbers without a 
single prime among them. And yet if we were to make the impossible- 
looking demand for a million consecutive composite numbers fac- 
torials provide the answer immediately. It is only necessary to 
write 1000001! + 2, + 3,. . . + 1000001, and we have two sets of 
numbers which fulfil the requirement. For since n! is divisible by 
2, 3, 4,. . ., adding or subtracting 2, 3, 4, etc., will not affect this 
property and all integers within the intervals proscribed must there- 
fore be composite. 

Of course this is a most profligate means of ensuring that such a 
demand is met. For instance to make sure of say, twenty-four con- 
secutive composite numbers it is necessary to start with 25! or 
15511210043330985984000000 whereas all the integers between 
3137 and 3163 are composite. 
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(3) A most interesting and far-reaching generalisation of Fermat’s 
Theorem was first announced by Euler in 1760. In this he began by 
examining those integers which are ‘relatively prime’ to one another, 
that is which have no common factors. Thus, for example, 7, 8, 9, 25, 
are all relatively prime though three of them are composite. 

Now any given number m has just so many integers less than m 
which are relatively prime to it; the number of such integers is 
denoted by ¢(7m)—Euler’s function. When m is a prime then clearly 
¢(m) = m — 1, but with m composite there is some counting to 
be done. Thus when m = 14, three are six integers prime to m, 
namely, 1, 3, 5,9, 11, 13, and therefore ¢(14) = 6. 

Euler’s generalisation in effect substitutes ¢(m) for (p — 1) in 
Fermat’s statement that n?-! — 1 is divisible by p, and extends 
the theorem to include all values of γι. We now have, “n®™) — | 
is exactly divisible by m, provided only that n is relatively prime 
to γ΄. 

Taking again the above example, m = 14; d(m) = 6, the theorem 
shows that for n = 1, 3, 5, 9, 11, 13; πθ — 1 = 14k. 


Thus: 39 — 1 = 728 = 14.52 
58 — 1 = 15624 = 14.1116 
96 — 1] = 531440 = 14.37960 etc. 


It should be noted that can be any integer that is relatively prime 
to m, and thus we have for instance, 


176 — 1 = 24137568 = 14.1724112. 


(Note. Since 7 is nota member of the set, 7° — 1 = 117648 + 14k.) 

Although quite incidental to the above ideas it will be of interest 
to observe that if any two or more of the integers 1, 3, 5, 9, 11, 13, 
are multiplied together and the product divided by 14, the remainder 
is always a member of the set. 


For example, 5.13 = 65 = (4.14) + 9, 
3.11 = 33 = (2.14) + 5, ete. 


This is a property of all sets of integers which are relatively prime 
to a given number. 

Methods of determining ¢(m) for any value of m can be stated 
simply. Their derivations and proofs will not, however, be given 
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here because they are most elegantly presented by the use of con- 
gruence techniques. (See H. Davenport’s The Higher Arithmetic.) 
The basic principle can be described as follows: 

Let m=a’.b*.ct...., where a, δ, and c, etc., are prime 
factors. 

Then ¢d(m) = m(1 — 1/a)(1 — 1/b)(1 — I/c). . .. (Note that r, 
5, t, etc. do not enter into the calculation.) Thus, for example, when 
m = 40 = 28 x 5, φ(40) = 40(1 — 1/2)(1 — 1/5) = 40 x 1/2 x 4/5 
= 16. 

In the following list (Table 13) all the values of ¢(m) for m < 400 
have been tabulated. Before looking at its practical use as an aid to 
factorisation the following observations can be made; all but the 
final one can be proved and can be used to extend the table as 
required. 


(a) Every value of ¢(m) is an even number. 

(δὴ) When m is prime, then ¢(@) = m — 1. 

(c) When m is odd, then ¢(m) = 4(2m). 

(4) When m is even, then 4(2m) = 2.9(m). 

(ὁ) When m is odd yet not a multiple of three, (i.e. of the form 
6k + 1), then $(3m) = 4(4m) = ¢(6m). 

(f) When m is a square number, say a?, then (a?) = αιφ(α). 

(g) All values of (7m) are repeated at least once. This is a con- 
jecture which has not yet been proved. 


The practical use to which this function can be put is perhaps best 
illustrated by taking an example from Table 13. It will be seen, for 
instance, that 6(399) = 216. (See the end of the last line in the table.) 
Since 10 is prime to 399, we can therefore state that 10716 — 1, 
(= 9.J216) is a multiple of 399, (= 3.7.19). Looking further up 
the table it will be seen that 216 also appears as a function of 
351, 333, 327, 259, and 247. Factorising these numbers we obtain 
the ‘new’ primes 13, 37, 109 and it is then possible to state that 
Igig is a multiple of 3.7.13.19.37.109. (The reader might find a 
gentle exercise in checking these findings back to the original 
definition of Euler’s function.) It will be noted that Fermat could 
not have helped here, since 217 is not prime. 


(4) It was recognised early in mathematical history that no formula 
could be found that excluded all composite numbers but there still 
remained conjectures about ‘forms’ of numbers which (a) provided 
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only primes, and (6) contained an infinity of primes. Perhaps the 
best known guess in class (4) was that made by Fermat. His formula 
for primes, to this day known as ‘Fermat Numbers’ or ἔμ was 
22" +1, where ἡ = 1, 2, 3,. . .. This expression clearly produced 
primes for n = 1, 2, 3,4, and there the matter stood until in 1732 
Euler showed that Fs = 4,294,967,297 = 641. 6,700,417. 

Nearly 150 years elapsed before Fg was cracked. 


Ἐς = 18,446,744,073,709,551,617 = 274,177 x 67,280,421,310,721. 


Many more composite numbers of this form have since been found 
but not, so far as I know, any more primes. Thus F, for n = 7, 8, 
10, 13, 14, 16, is known to be composite, and the following factors 
have been found: 


Fy is divisible by 37.216 + 1 

F\, is divisible by (39.218 + 1)(119.218 + 1) 
Ἐπ is divisible by (7.214 + 1)(397.216 + 1) 
Fis is divisible by 579.271 + 1 

Fis is divisible by 13.229 + 1 

F 23 is divisible by 5.275 + 1 

F3¢ is divisible by 5.289 + 1 

Fg is divisible by 3.241 + 1 

F73 is divisible by 5.275 + 1 


(To appreciate the size of these numbers the last factor, 5.275 + 1 
is equal to 188,894,659,314,785,808,547,841.) 

These numbers increase at a fabulous rate and yet, thanks to 
highly sophisticated techniques developed in the last half century, 
it is now known that the almost incomprehensible number Fig45 is 
divisible by 5.21947 + 1. 

It would seem that Fermat’s conjecture has somersaulted and one 
should now ask whether any primes beyond F4 exist at all. 


(5) Although not strictly belonging to either of the above classes, 
no discussion of the primes would be complete without some 
mention of the ‘Mersenne’ numbers. Mersenne, a contemporary of 
Fermat (1588-1648), studied the numbers 2? — 1 (p—prime) and 
in 1644 made a statement to the effect that the only values of p 
- that make 2? — 1 ἃ primeare 2, 3, 5, 7, 13, 17, 19, 31, 67, 127, 257. 
(Being nobody’s fool one assumes that he added, or implied, 
p < 260.) 
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These numbers do not increase with anything like the rapidity 
of Fermat’s, but nevertheless the factorisation of the intervening 
numbers must have presented, in those days, a most formidable 
problem. It was not until the end of the nineteenth century that 
Mei was proved to be prime and this was followed by Cole’s 
demonstration in 1903 that Me7 was composite (193,707,721 x 
761,838,257,287). More primes of this form, p > 257, have been 
found since, the largest to date (1968) being 211715 — 1. 

Mathematical fortune-telling on meagre empirical evidence is no 
longer a popular pastime but it must not be forgotten that the efforts 
to prove or disprove these two wild guesses have probably contri- 
buted more to the development of factorisation methods and tests 
of primality than any other problems in the Theory of Numbers. 

In class (6) we have a rather different situation. Since, for instance, 
4n +1 and 6” +1 contain respectively all the odd numbers and 
all those not divisible by three, these progressions obviously contain 
all the primes. Simple as it may seem the proof that each of these four 
progressions contains an infinity of primes requires some effort 
and any amateur who completes the set can call himself a mathe- 
matician. 

Dirichlet proved (1837) that the series an+b5 (a relatively 
prime to 5) generates an infinite number of primes and it is con- 
jectured—as yet without proof—that the same applies to the series 
n? +1. There is, of course, no limit to propositions of this kind 
but when it comes to proof it is difficult to find a point of contact 
with the few known properties of the primes. 


EXERCISES 
1. For what values of n is n® — 1 divisible by 20? 


2. Use Table 13 to find four prime factors of 512 — 1. 
3. What is the value of φ(402) 


Ἀ PRIME NUMBERS (Part 2) 


(1) In this chapter some of the lesser known surmises and observa- 
tions concerning the primes will be described. 

One of the most curious of these is Euler’s proposed test for the 
primality of integers of the form N = 4n + 1 whichendin 3 or 7, 
or more concisely those of the form 20m - 13,+ 17. (These 
integers represent one quarter of all the odd numbers.) The criterion 
then runs: 

“Let R be the remainder after subtracting from 2N the next 
smaller square which ends in 5, namely (5n)?. Then to R add 
cumulatively the numbers 100(m — 1), 100(m — 3), 100(m — 5),.. . 
and so on. If among these sums, including R itself, there occurs a 
single square number, then N is either a prime or is divisible by 
this square. If, on the other hand, two or more squares, or none, 
occur then N is composite.’ 

The following four examples will serve to disentangle these unusual 
and somewhat confusing instructions. 


(A) N = 637. 2N = 1274 
352 = 1225 = (5n)?, son = 7. 
R= 49 


Then R, (R + 600), (R + 400), (R + 200) = 49, 649, 1049, 1249. 
In this set the only square is 49, and therefore N is either prime or 
is divisible by 49. In fact 637 = 49-13. 


(B) N = 2437. 2N = 4874 
652 = 4225 therefore n = 13 
R= 649 


The sequence then becomes, 649, 1849, 2849, 3649, 4249, 4649, 4849. 
Here the only square number is 1849 and since 2437 is clearly not 
divisible by 1849 then it is prime. 


i 
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(C) N = 2303. 2N = 4606 
652 = 4225 n=13 
R= 381 


Sequence—381, 1581, 2581, 3381, 3981, 4381, 4581. Since this does 
not contain a square, N is composite. 


(D) N=2117. 2N = 4234 
65? = 4225 n= 13 
R= 9 


Sequence—9, 1209, 2209, 3009, 3609, 4009, 4209. Of these numbers 
9 and 2209 are square and therefore 2117 is composite (= 29-73). 

One might go on for a long time testing numbers of these two 
forms at random before coming across an anomaly. In fact where 
the test indicates that an integer is composite it is invariably true 
but unfortunately some composite numbers appear amongst those 
it picks as prime. The first few of these are 153, 333, 477, 657, 
833,. .., and are in general those integers which are multiples 
of primes of one of the forms (20k + 13) and (20k + 17) and 
also of one or more of the squares of 3, 7, 9, 11, 23, 27,. . . @e. 
(20m + 3, 7, 9, 11)?). 

(Assuming that the first step in any test for primes would be the 
elimination of 3,7, 11, 13, as possible factors the first error would 
not occur until we reached 8993 = 17.23%, a truly remarkable 
achievement.) 

Euler’s test for primes of the above forms can therefore be polished 
by the further provision that any integer appearing as a prime must 
again be tested for divisibility by the numbers (20m + 3, 7, 9, 11)?. 
In passing it might be of interest to note that the numbers 1353. 17, 
132 . 37, 132. 57,. . . have sequences which include three squares, 
whilst 32. 13, 72. 13, 92. 13, 112. 13,. . . provide the cases where 
there is one square which divides N. 

As far as I know no one has put forward anything comparable 
as a test for the other three forms of odd numbers and this most 
peculiar approach to the theory of primes—with such a master as 
its sponsor—might well initiate a rewarding project for some amateur. 


(2) The following miscellaneous criteria for primes are given without 

references; they have been gathered much as one picks blackberries 

and I cannot say with certainty that all of them have been proved. 
6 
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They should, however, give some idea of the concentrated effort 
and the diversity of approach which has been directed upon these 
fascinating and elusive members of our number system. Perhaps some 
of them will provoke a renewed attack. In what follows the symbol 
(if) will be used to indicate the expression ‘if, and only if”. 


(a) Anumber is prime (if) it is not expressible in the form ab + xy 
when a,b,x,y, are positive integers such that a+ b= x — y. 

(ὁ) A number p is prime (if) it occurs (p— 1) times in the 
(p — 1)th set where the first set is 1, 2, 1; the second set is 
1, 3, 2, 3, 1—formed by inserting between each two terms of 
the preceding set their sum—the third 1, 4, 3, 5, 2, 5, 3, 4, 1, 
and so on. 

(c) If πὶ is an odd number then 4n + 1 is prime (if) it is a factor 
of (22 + 1)/5. 

(d) If p is a prime of the form 4k+1 then 2p +1 is also a 
prime (if) it exactly divides (2? + 1)/3. 

(e) If p is a prime of the form ἀκ —1 then 2p +1 is also a 
prime (if) it exactly divides 2? — 1. 

(7) If p is an odd prime then the integral part of (4/5 + 2)? — 
2?+1 is exactly divisible by 20p. 

(5) A number of the form 6n+1 is prime (if) m cannot be 
expressed in either of the forms 6xy + x + y, or 6xy — x — 


7. 

(h) Anumber of the form 6n — 1 is prime (if) 7 is not of the form 
6xy +x—/y. 

(i) If p and gq are different primes then ρας} + q?-1 —1 is 
exactly divisible by pq. 

(j) At least four primes lie between the squares of two consecutive 
primes (>3). 

(k) If 4n —1 and 8 --- are both primes, then 24-1 —1 isa 
multiple of 87 — 1. 

(ἢ Gf) p is prime then 141/24+1/3+...1/—~—1) is 
exactly divisible by p. 

(m) Given that s and ὃ are respectively the arithmetic means of 
the squares and cubes of all the numbers less than and rela- 
tively prime to n, then n? = 6ns — 4c. 

(n) Between any two primes there is an odd number of composite 
integers. 
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(3) Although nobody nowadays searches for a formula which will 
generate nothing but prime numbers it is still reasonable to ask 
‘what is the number of primes (say P(x)) less than a given integer 
x?’. Most surprisingly one answer was given in 1896 by Hadamard 
and de la Vallée Poussin independently, both men in or about 
their thirtieth year. It was to the effect that P(x) approaches the 
value x/logex as x tends to infinity. 

This intrusion of Analysis—which deals with continuous functions 
and might be called the science of approximation—into a branch of 
mathematics concerned exclusively with the whole numbers was 
almost as revolutionary as the introduction of the minus sign into 
primitive arithmetic. 

Of course for some limited values of x, P(x) can be determined by 
an actual count, and the following lists may be worth noting for 
reference. 


Number of primes 


0 --, 1000 168 
1000 — 2000 135 
2000 — 3000 127 
3000 — 4000 120 
4000 — 5000 119 
5000 — 6000 114 
6000 — 7000 117 
7000 — 8000 107 
8000 — 9000 110 
9000 — 10000 112 


On a more extended scale we have (excluding 1, and 2) the 
number of primes less than 


10 = 3 

100 = 24 
1000 = 168 
10000 = 1229 
100000 = 9592 


10000000 = 664579 
100000000 = 5761455 
1000000000 = 50847478 
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(4) Finally a word on large primes. About sixteen years ago 2137 — 1, 
a number of some forty digits, held the record. Shortly afterwards 
180(2127 — 1)? + 1 (about 83 digits) was proved prime. Since then 
we have had 2? — 1 where p = 521, 607, 1279, 2203, and 2281, 
this last containing 686 digits. 

For a time it was thought by some that 2? — 1 was prime when q 
is a Mersenne prime but it has recently been shown that 22-1 — 1 
(= 28191 — 1) is composite. And now there is this fabulous prime 
211213 — ], employing 3376 digits. 

Even this may have been surpassed by the time these notes appear. 


AN INTRODUCTION TO 
DIOPHANTINE EQUATIONS 


(1) Equations involving two or more unknowns which require 
solutions in integers or rational fractions date back to Diophantus of 
Alexandria in about the third or fourth century A.D. It has taken 
almost all the intervening time—some fifteen centuries—to bring 
these problems under control of some sort and even now there still 
remains the question of whether x” + y" =z", for n greater 
than 2 has any solution. 

Linear Diophantine equations can be expressed in the form 
ax + by = cc, in which a, δ, and c¢ are given and it is required to 
find integral values of x and y. 

We are all familiar with the type of popular mathematical puzzle 
exemplified by the following: A farmer pays £7 each for sheep and 
£12 each for calves. (I wouldn’t know whether these are realistic 
prices today.) If he buys some of each for a total bill of £71, how 
many sheep and calves has he acquired? 

Taking the number of sheep and calves to be respectively x and y, 
we then have the equation 7x + 12y = 71. It is to be understood 
of course that at the end of the transaction the animals are to be 
alive and running about on all four legs. 

The standard school-book treatment of such a problem proceeds 
as follows: 


7x + 12y =71 (1) 
Dividing by the smaller coefficient (7) we have, 
x + y + Sy/7 = 10 + 1/7 
or x+y + Gy — 1/7 = 10 (2) 


Since x and y are integers then (Sy — 1)/7 must be an integer 
and consequently (15y — 3)/7 is also an integer. (Here a multiplier 


| & 
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is chosen which will make the coefficient of y differ from a multiple 
of 7 by unity.) 


Now (15y — 3)/7 = 2y + y/7 — 3/7 


= 2y + (y — 3)/7 
As before it follows that (Ὁ — 3)/7 is an integer and equal, say 
to ἡ. Then y—3=7n or y=7n+3 (3) 


Substituting this in (2) we have: 


x + In +3 + (35m + 15 — 1/7 => x+ 12n+5=10 
or x = 5 — 12a. 


For x to have a positive value ἡ must be zero and we have x = 5, 
and from (2) y =3. 


(2) Now there is a much simpler method of arriving at the above 
result. The algorithm may take a little space to describe but, once 
understood it is much quicker and more reliable in dealing with 
equations having large coefficients of x and y. 

Taking again the above example we re-write the equation thus: 
7x = 71 — 12y. (This is not really necessary but it helps with the 
explanation.) Divide the coefficients on the right by 7 and note 
down the remainders, i.e. 1(c) and —5(b). (Remember the general 
equation is ax + by =c.) Set down the numbers 1, 2, 3,... 
(a — 1), or in this case 


Τὰ 9, My OO 


Multiply each of these by the ‘b’ coefficient, namely 5, to give 
5, 10, 15, 20, 25, 30. Divide each of these by 7 and list the 
remainders thus: 


3, 3». 4, 64,2 


Note that the ‘c’ coefficient, namely 1, is at the third position and 
write y = 3. (As easy as that) x = 5 is easily found by substitution. 

The process is far simpler to apply than to describe and as a further 
example of the power of this method we ask ‘What are the smallest 
integers which will satisfy the equation 


l1x + 127y = 1,067,723? 
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Rewriting we have 11x = 1067723 — 127y. Dividing by 11 the 
remainders are 8, 6. Then 


1, 2, 3, 4, 5, 6, 7, 8, 9, 10. 
And, multiplying by 6 and subtracting multiples of 11— 
G54, 1, δειναῖς 
Since 8 is at the fifth position then y= 5, and by substitution 
x = 97008. 


Check: 11 x 97008 = 1067088 
xX hy τῷ 635 


1067723 


(3) For equations of higher degree the problem becomes more 
difficult and although some types are amenable to treatment there is 
in fact as yet no general technique for determining even whether a 
solution is possible or not. 

In the case of equations describing the so called Pythagorean 
right angled triangles we are on familiar and well trodden ground. 
These are triangles of sides x, y, and (hypoteneuse) z, and the 
problem is to find solutions of the equation x? + γὃ = ζῇ, x, y, z, 
being integers and having no factors in common. 

The general solution can be reached by the following reasoning— 
I give this in full because it indicates the best known approach to the 
more difficult problems in this terrain: 


(1) x and y cannot both be even, for then z would also be even. 
(2) x and y cannot both be odd for then x? + y? would be of the 
form 4n + 2 which is never a square. 


Suppose then that x is odd and y, even. We can now write the 
equation x? + 413 = z? or, 413 = (z + x)(z — ). 

But x and z are odd and therefore z+ x, and z—~-x are 
both even. Putting z+x=2s and z—x=2r we then have 
413 = 4sr, or 13 = rs. 

Also (z + x) + (2 — x) = 2s+2r or z=r-+s and similarly 
x=s—r. 

Since x and z are relatively prime then so are r and 5. In 
consequence the equation 13 = rs requires that r and s are each 
perfect squares, say r=? and s=m?, Then t= mn and it 
follows that 

x =m? — n?; y = 2mn; 2 = m? + ηϑ. 
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Even so m and nm must obey certain conditions. Thus m is 
greater than m so that x is positive; m and nm must have no 
common factor or it would be shared by x, y, z; m and n must 
not both be odd for then x and z would both be even. 

With these reservations in mind it is possible to tabulate as many 
solutions—in integers—of the equation x? + y? = z? as we wish. 

For instance the first few are: 


m 


AnNUUhARWH 
-—UANReEWNHeE B 
-- 
S8BSaoRE5a & 
_ 

Δ 


5 42 3] 


The above general solution to the Pythagorean equation has been 
well documented but there is a geometrical consequence to this 
which is perhaps not so well known. Given a right angled triangle 
whose sides and hypoteneuse can be expressed in integers we consider 
the radius of the inscribed circle. 


The area of AABC = x/2y = As AOB + BOC + COA. There- 
fore x/2y = rz + dry + rx = (x + y +2) and r = xy/(x+ 
z). 
᾿ πον we replace x, y, Ζ by their equivalents in m and n we 
have 
(m2 — n®) . 2mn 
= n(m — n) 


P= 


-- ...........-.-. 
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Since both m and n areintegers then r, the radius of the inscribed 
circle is also always an integer. In particular when x = 3, y = 4, 
z= 5, then r= Il. 


EXERCISES 


1. Solve 5x + 9y = 1001, for x and y. 


2. If in the equation x? + y? =z? we give z the value of 26, 
what are the corresponding values of x and y? 


10 SOME OBSERVATIONS ON Χ' 


(1) In a broad sense experiments are of two kinds. Those which are 
made with the object of confirming or disproving a theoretical idea, 
and those whose sole purpose is that of gathering information. 
Of the two the latter can usually be expected to provide more oppor- 
tunities of excitement than the former (unless of course the theory 
happens to be particularly novel and one’s own). 

Experiments designed to tell us how numbers behave under certain 
specific conditions nearly always produce something of interest 
if they are carried far enough. And in the world of numbers it is 
indeed remarkable how little initial data is needed to start a train of 
successive observations. 

Let us take as a starting point, for example, one of the simplest 
properties of integers in the scale of ten. ‘All powers of 10 end with 
the “units” digit 0.’ (Incidentally all powers of x inthe scale of x 
do the same but we shall confine ourselves to the everyday system 
here.) I am sure that no reader will require a formal proof of this. 
Nor will he need telling that all powers of 5, or for that matter 
of any number of the form 10k + 5, end in 5 and that those of 
6 end with a 6. The same thing holds for numbers terminating 
in 1 and it is easily seen that the even powers of 10k + 9 endin 1 
and the odd powers in 9. 

Having made these observations it is natural to ask next what 
happens to the ‘units’ digits of powers of the remaining numbers 
10k + 2,3,4,7, and 8. Since we are only concerned at the moment 
with the final digits it is a simple matter to construct a table without 
calculating out the powers in full. For instance the units digits of the 
powers of 2 are 2, 4, 8, 6, 2, 4, etc., those of 3 are 3, 9, 7, 1, 3, 
and so on. In fact it is quickly seen that the terminal digits of the 
powers of all the natural numbers repeat themselves after, at most, a 
cycle of four digits. 
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Table 14a 
Terminal digits of x” 
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0 
1 
2 
3 
4 
5 
6 
7 
8 
9 


Having produced this table, what can be learnt from it? For a start 
there are four immediate observations which can be made. 


(A) All the odd powers of any number x may end in any of the 
digits 0, 1,2, 3,. . «9. 


(B) x, x8, and in general x4”*? can only end with the digits 
0, 1, 4, 5, 6, and 9. 


(C) x is restricted to endings of 0, 1, 5, 6. 


(D) The even powers of any number can never terminate with the 
digits 2,3, 7, or 8. 


Much more information, however, can be squeezed out of this 
table. From at least the time of Fermat and Mersenne, numbers of the 
general form x®+1 have been the object of much study and 
speculation particularly in their connection with prime numbers. 
Although the table has no direct bearing on the primes it does enable 
us to sweep a wide area free from one class of composite numbers. 

Since all integers ending in 0 or 5 are divisible by 5 then 
clearly the addition of unity to those powers of x which have endings 
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of 4 or 9 will make the resulting numbers composite (i.e. multiples 
of 5). Similarly, subtracting one from powers ending in 1 or 6 will 
have the same result. 


Thus 2? + 1, 26+ 1,. . .24*% + 1, 

and 25 — 1,28 —1,. . . 2% — 1, are all divisible by 5. 

Multiples of five are of course easily discernible when expressed 
in integers but this simple table enables us to list all such multiples 


which are of the general form x” +1 in the following specific 
algebraic statements: 


(10m + 2)4%*—1, (10m + 2)4*+2 + 1, 
(10m + 3)4%*—1, (10m + 3)4*+2 + 1, 
(10m - 4)**—1, (10m + 4)?#+1 + 1, 
(10m +6)" —1, (10m + 7)**+2 + 1, 
(10m + 7)4*—1, (10m + 8)*+2 + 1, 
(10m + 8)4*—1, (10m + 9)4*+1 + 1, 
(10m + 9)4* — 1. 


(2) Now in considering the terminal (i.e. units) digits of x" we 
have, in effect, been concerned only with the remainders left after 
division by ten. It will be reasonable to continue with a similar 
examination of the ‘digital roots’ of the powers of the natural num- 
bers; these are of course the remainders left after division by nine. 

The corresponding table is quite simple to construct; as before it 
is not necessary to calculate the actual powers of x since if x = 
9r+ a, then x? = 9s + a®, χϑ = Θὲ + a®, etc. Consequently we 
need only multiply the digital root of x successively by a. 

For example when x = 7, the digital roots of x, x2, x3,... 
are found thus: 7, 7:7 = 49 (D.R. = 4), 74=28 (D.R. = 1) 
and so on. In this manner we arrive at Table 14b which again is 
quite general and repeats itself cyclically when extrapolated either 
vertically or horizontally. 

Since all numbers having the digital roots 3,6, and 9 are divisible 
by three, we can use this table to sieve out another set of composite 
values of x” + 1. 

We see, for example that the digital roots of 21, 23, 25, 27, etc., 
are always one of the digits 2, 5, 8, and that those of 23, 24, 26, 28, 
etc., are either 1, 4, or 7. It follows therefore that 3 is always a 
factor of both 

22k+1 + 1 and 22% — 1, 


ws 
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Table 14b 
Digital Roots of x” 
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1 
2 
3 
4 
5 
6 
7 
8 
9 


Proceeding on these lines and generalising as before we are able 
then to make a list of all the values of x” + 1 which are multiples 
of 3. Thus: 


(9m + 2)?*—1, (9m + 2)**+1 + 1, 
(9m - 4} —1, (9m + 5)**1 +1, 
(9m + 5)*—1, (9m + 8)*1 +1. 
(9m - 7) —1, 
(9m + 8)** — 1. 


(3) These two tables, extremely simple both in construction and 
appearance, contain a surprising amount of information, some of 
which is not easy to uncover by ordinary algebraic processes. To 
take an elementary example let us assume we are searching for primes 
among numbers of the form 2” + 1. Table 14b tells us that 274+1 + 
1, an expression containing all the odd powers of 2, is divisible by 
3, whereas it is seen from Table 14a that 24*+2 + 1, or 22-6-10,ete 4 
1, is always a multiple of 5. Thus we can only expect to discover 
primes among numbers of this form when n = 4k, i.e, 2% + 1. 
This news will not startle Number Theory addicts but it is still a 
fair stride from our initial observation that ‘all powers of ten have 
the terminal digit 0’. 


ao ii .“ὦὃὦὃὕὕ 
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The enquiry started at the beginning of this chapter can, of 
course be carried much further and indeed the examination of the 
remainders left after dividing x” by the successive primes 7, 11, 
13,. . ., leads into an extremely interesting and advanced field. 

However our two primitive tables (14a and 146) can still provide a 
mass of information about exponential numbers which are multiples 
of 5 and 3. 

It will be seen at once, for example, that x” + x"*? always has 
the terminal digit 5 (and hence is a multiple of 5) when x = 10k + 
2, 3, 5, 7, 8. Similarly 25 + 2511 is obviously a multiple of 3 and 
this can be extended to x” + χη when x = 9k + 2, 5, 8. 

The more general numbers x” + y” can also be separated into 
multiples of these two factors and a few examples are given here which 
the reader may care to disentangle for himself. 


Multiples of 5. 


Q4ar+l 4. 348+1 
Q4ar+l 4. 748-1 
Q4r-1 4. g4s—1 


and as an example of still wider generalisation, 
(10k + 7)#*1 — (107 + 3)4-!. 
where k, j, r,s, take any of the values 0, 1,2,3,.. .. 


Multiples of 3. 


Qor+1 4. 568+2 
qntl — Jn 


and (9k + 2)8r-1 + (9) 4 5)ss+2, 

The use of these tables is not, of course, confined to the addition 
or subtraction of the powers of only two integers. Thus, for instance, 
if we consider the horizontal line of digits against x = 2, in Table 
14a, we have, 


2+44+8+6+2+4= 26, 

and it follows that 
(2 + 22 + 25 + 24 + 25 + 2%) — 1 

is divisible by five. 
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In a similar manner Table 14 tells us that 
19 4+ 21 + 32 + 43 4 54 + 65 4 76 


is a multiple of nine. 

The above principles can be applied in a variety of further instances, 
not all of them trivial, but the point I am hoping to have made is that 
sometimes the most unexpected findings can emerge from the develop- 
ment of the simplest of elementary observations. 


EXERCISES 
1. What is (a) the ‘unit’s’ digit, and (6) the digital root of 1111} 
2. In the progression 3°, 31, 32, 33,. . . 3", what are the powers 
which have the remainder 1 after dividing by 7? 


3. Given that the sum 18 + 23 + 33 + 48 + x3 is a multiple of 
nine, x (<10) can have two values. What are they? 


‘ f FINITE ARITHMETIC 


(1) In the last chapter we examined a particular class of numbers— 
the powers of x—with special regard to their remainders after 
dividing by both 9 and 10. Indeed we started from two simple 
truths which were obvious by definition, namely that (1) the units 
digit of any integer remains the same whatever multiple of 10 
is added to it, and (2) the ‘digital root’ of any integer is constant 
for all additions of 9k. 

We dealt there with two special cases of what is clearly a general 
statement. If we are considering only the units digits of the integers 
it might be said that 0 is synonymous with 10, 20, 30, etc., or 
among digital roots 0 is identical with 9. Similarly it would be 
just as convenient to mark Ὁ on a clock face instead of 12, since 
the cycle starts again at this point. Thus the clock will register the 
same time in twelve hours or, with perhaps some astronomical or 
horological reservations, in a thousand years, or any numbers of 
revolutions. Or to give a further example, whenever two values of n 
differ by an even number then (—1)" remains identically the same. 

This notion of the ‘sameness’ of numbers which differ only among 
themselves by some specified multiple was first crystallised and 
given a formal mathematical symbol by Gauss towards the end of the 
eighteenth century. Incidentally, he was in his teens at the time. 

What Gauss’s notation expresses in fact is that if two integers, 
a and ῥ, differ by a multiple of a particular number m then a is 
said to be congruent to 6 with respect to the modulus m. In this 
notation we write ὦ ξξι (mod m) instead of the more familiar 
algebraic a—b=xm. Arithmetical congruence then implies 
equality except for the addition or subtraction of some multiple 
of m. It will immediately be clear that in the arithmetic ‘modulo 
m’ there are no integers greater than m—1, hence the term 
‘finite arithmetic’. 

Congruences have much in common with ordinary equations 
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and indeed are manipulated according to exactly the same rules 
with the one proviso that the same modulus is retained throughout 
the operations. 

Thus, if a = x (mod m), then, all to (mod m), we have 


ar = xr 
at+r=x-4+r | 
a’=*x' etc. 
If in addition b=y 
then ab = xy 
and at+b=x+y 


Although a congruence can always be multiplied throughout by an 
integer, cancellation of a factor is not permissible unless the factor is 
relatively prime to the modulus. For instance it is legitimate to 
reduce the congruence 48 = 18 (mod 10) by a factor of 3 to 
give 16 =6 (mod 10), but we cannot divide by the even factor 6 
(which is not prime to 10) without arriving at the false result 
8 =3 (mod 10). 


(2) The concept of congruence is no mere mathematical toy; 
it enables us to express many propositions with economy and ele- 
gance and particularly to operate with numbers so large that ordinary 
methods of calculation would be impossible to apply. 

As an elementary example of this technique it might be interesting 
to re-examine the ‘extraction’ of digital roots. Since any integer N 
can be written do + 10a; + 102a2 + 10843 +... + 10"a,, and 
10 =1 (mod 9)—from which it follows that 10?=1, 108 =1, 
...10%=1—it is clear that N=ap+ai+az2+as+... 
+ a, (mod 9). Hence N is equal to a multiple of 9 plus the sum 
of its digits and if the latter is divisible by 9 then so is N. 

Similar reasoning is equally effective in proving the common test 
for divisibility by eleven. 

For the benefit of any newcomers who find the technique 
unfamiliar perhaps the following examples will help to clarify the 
primary methods of operation. 


(A) We have 9=2 (mod 7), 


then 92=22=4 
7 
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and 98 = 36 = 1, and therefore 9% divided by 7 
leaves the remainder 1. 


(B) Since = —1 (mod 7) 
62 =(-—1)? =1 

and hence 62" == 1 (mod 7) 

(C) (mod 7) 


and since 15=1 (mod 7), and 36 =1 (mod 7), it follows that 


5 = 1/3 (mod 7) 
3 = 1/5 (mod 7) 
6 = 1/6 (mod 7) 


It is clear that this arithmetic, although confined by the modulus, 
is free to operate with both negative signs and fractions. 

As an example of the power of the congruence technique we will 
now apply it to an actual computation. It is required, let us say, to 
find the remainder upon dividing 539 by 103. Remembering froma 
previous chapter that any number can be expressed as the sum of 
single powers of 2, we first note that 526 = 52 x 58 x 516, 


With this in mind we now proceed as follows: 


5? = 25 (mod 103) 


54 = 625 =7 
58 = 49 
516 = 2401 = 32 


It is now required to calculate 53 x 58 x 510 (mod 103) and from 
the above figures we have, 


55 x 58 = 25 x 49 = 1225 = 92 
92 x 516 = 92 x 32 = 2944 = 60. 


60 then, is the remainder after dividing 535 by 103. To reassure 
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those who might doubt so facile a result, here is the actual division 
done the hard way. 


103/1490116119384765625 
103 


460 481 446 
412 412 412 
481 699 345 
412 618 309 
691 813 366 
618 721 309 
736 928 572 
T2L S27 818 
151. 14} 575 
103 103 515 


481 446 60 = remainder. 


(3) It will be of interest at this point to demonstrate the strength of 
modulus arithmetic by an example of its application to a principle 
of far reaching importance in Number Theory. Fermat’s theorem 
states that if a is any integer and p a prime, then a? —a is 
divisible by p. Or in congruence notation that a? =a (mod p). 
Nearly twenty five centuries ago the Chinese knew this to be true in 
the case of a = 2, and until the early nineteenth century the converse, 
namely that such a factor p is essentially prime, was accepted. 
Then in 1819 Sarrus discovered that 2341 — 2 is divisible by 341, 
a composite number (11 x 31). 

The actual computation of the integer 241 — 2, containing over a 
hundred digits, and the subsequent tests for divisibility would present 
a formidable task, but see what light work our modulo arithmetic 
makes of it. 


We have first 
Then it follows that 


25 = —1 (mod 11) 
(25)2n = 1 


and hence (25)68 = 2340 = ] 
Similarly 25 = 1 (mod 31) 
and 2340 = | 


98 Finite Arithmetic 


Thus both 11 and 31 are divisors of 2340 — 1 and hence also 
of 2341 — 2. 

Fermat’s theorem is perhaps more usually expressed in the terms, 
‘if p is a prime and a is not a multiple of p then a?-!—1 is 
divisible by p’. (More refined statements of the theorem will be 
referred to later.) 

That the converse of this is not always true can be seen from a 
much simpler example than the above. Putting a = 4, and p = 15, 
(a composite number) we have 444 — 1 =0 (mod 15). 


Now = —1 (mod 5) and 4=1 (mod 3) 
So 414 =] 414 =] 


Hence 444— 1=0 (mod 5) and 434 -- Ἰ Ξε 0 (mod 3). Asa 
check we have 414 — 1 = 268435455 which is obviously divisible 
by both 3 and 5. 

There are other variants of this famous theorem which lead into 
much more advanced theory than can be embarked upon here. 
As examples I will merely comment that the two following statements 
are true: 


3120 — 1 = (0 (mod 112) 
21092 _ 1 =0 (mod 10932). 


(4) Turning now to another development of congruence theory 
we consider the squares of the natural numbers, 1?, 25, 3%,.. .. 
Listing the ‘digital roots’ of these successive squares we have: 


n® 149 623 56 @ δὲ 81 WO PI. :ς 
mee he OO OF Pi BA a ΨΥ a uae 


The sequence is clearly repetitive and indeed it is easily shown 
that however far it is continued no other digits than 1, 4, 7, 9, 
appear, and since in fact digital roots represent the remainders left 
after dividing by 9 the set might just as well be written 0, 1, 4, 7. 
These digits are known in congruence language as ‘quadratic residues’ 
modulo 9. 

It is easily shown that a closed cycle of this sort will occur with 
any divisor. Let the remainders (residues) on dividing the successive 
squares by πὶ be ri, re,rs,.... Thenif x is any number we have 
(nm + x)? = n? + 2nx + x? = x? (mod n) and ra +x will have 
the same value as 7 with a periodic recurrence of the terms. Using 
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a similar argument (m — x)? = x® (mod n) and hence ra — x = fz. 
Therefore for the modulus m there are not more than 4n, orif n 
is odd, 4(m — 1), residues to the square numbers. 

Quadratic residues to any modulus can be calculated quickly 
without actually writing down the sequence of square numbers. 
Since (x + 1)? = x2 +2x+1 it follows that ,+1=“%+ 
2x + 1 (mod n) and for example the residues modulo 13 can be 
derived in either of the following ways: 


>  ἊἮΝΖ rn x ἴῃ 


1 1 1 
2 2 4 
3 3 9 
4 4 (94+23+1)-13= 3 
5 5 (3424+1)-13=12 
6 δι: δ᾽ 10 6 (12+25+4+1)—13=10 
7 7 (10+2.6+1)—13=10 
8 8 (10+2.7+1)—13=12 
9 9 (12+28+41)—26= 3 
0 0 (34+29+4+1)—13= 9 


etc. etc. 


Integers which are not residues (mod m) are known as non-residue 
and the two classes are connected in the following manner. For any 
modulus n the product of two or more residues (r?, γα) is also a 
residue. (For if r? = x? and r4= y? (mod 7) then r?r¢@= (xy)? 
(mod n).) Also the product of two non-residues is again a residue, 
but the product of a residue and a non-residue is a non-residue. 

In the following table (Table 15) quadratic residues of the odd 
numbers 3,5, 7,. . .49, are listed. In order to show the distinction 
between a prime and composite modulus the table is divided into 
two sections. (The residues are given in the order of their appearance 
in the progressive sequence of square numbers.) 
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Table 15 
Quadratic Residues (to n?) 
n residues 
a 
, «aoe 
, bas 
ΕΓ 61,4,3,35,3 
iS ἘΦ 3, 2, 10 
17 = 1, 4, 9, 16, 8, 2, 15, 13 
19 1,4, 9, 16, 6, 17, 11, 7, 3 
2. 1,4,9. 16, 2, 13, 3, 18, 12, 8, 6 
29 1.4.9. 16, 25, 7, 2, 6,25, 13, 5, 28,24, 22 
31 +1, 4, 9, 16, 25, 5, 18, 2, 19, 7, 28, 20, 14, 10, 8 
3) «1, 4, 9, 16, 25, 36, 12, 27, 7, 26, 10, 33, 21, 11, 3, 34, 30, 28 
41 1,4, 9, 16, 25, 36, 8, 23, 40, 18, 39, 21, 5, 32, 20, 10, 2, 37, 
33, 31 
43 1,4, 9, 16, 25, 36, 6, 21, 38, 14, 35, 15, 40, 24, 10, 41, 31, 


29. 1,43, 11] 
47 1,44,9, 16, 25, 36, 2, 17, 34, 6, 27, 3, 28, 8, 37, 21, 7, 42, 32, 
24, 18, 14, 12 


9 1,4,0,7 
iS 189,14, © 6,4 
21 = 1, 4, 9, 16, 4, 15, 7, 1, 18, 16 
25 1,4, 9, 16, 0, 11, 24, 14, 6, 0, 21, 19 
“ 4, & δ ΣΎ oe 10, 0, 19, 13, 9, 7 
o> 1 9 10, Ὁ 3, 16, 31, 15. ἘΡΩ͂ ἘΣ, 4, Sh, 27,25 
35s 1, 4, 9, 16, 25, 1, 14, 29, 11, 30, 16, 4, 29, 21, 15, 11, 9 
39-1, 4, 9, 16, 25, 36, 10, 25, 3, 22, 4, 27, 13, 1, 30, 22, 16, 12, 10 
45 1, 4, 9, 16, 25, 36, 4, 19, 36, 10, 31, 9, 34, 16, 0, 31, 19, 9, 1, 
40, 36, 34 
49 1, 4, 9, 16, 25, 36, 0, 15, 32, 2, 23, 46, 22, 0, 29, 11, 44, 30, 


18, 8, 0, 43, 39, 37 


It will be seen from the above table that when 7 is a prime 
there are exactly 4(m— 1) residues, all of which are different 
integers. When n is composite the complete cycle still contains 
}(n — 1) residues but some of these are repeated one or more times. 
When n is a square or contains a square factor, 0 of course appears 
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and may be repeated several times in the complete cycle of residues. 
The number of different integers comprising the quadratic residues to 
a given modulus z, tell us then whether n is prime or composite. 
At first sight this may not seem a very practical method of establish- 
ing the primality of » but in fact it provides a very useful technique 
for detecting not only whether n is composite, but if so what are 
its factors. 

As a simple example let us take » = 13 x 37 = 481. Since all 
squares less than 481 are residues of n, we have Q.R.(m) = 1, 4, 
9,. . .441 and the sequence can be continued, since 441 = 213, 
and (2 x 21)+1= 43, 


441+43—481= 3=22? (mod 481) 
3 + 45 = 48 = 23? etc. 
48 + 47 = 95 
95 + 49 = 144 Stop. 


144 = 12? which has appeared before, and therefore n is 
composite. Further, since 144 is a residue of the 25th square 
we now have 252 = 12? (mod 481) or 252 — 122=0. That is 
(25 — 12)(25 + 12) = 13 x 37, the factors of 481. This will be 
dealt with more fully in a later chapter on factorisation methods. 


(5) The study of quadratic residues, to which the above is but the 
briefest of introductions, has produced many theorems of great 
importance in the realm of Number Theory but it would not be 
practicable to advance beyond this stage in a book of this sort. 

There are however some observations to be made on the residues 
to higher powers (than squares) which are neither difficult to grasp 
nor without their peculiar interest. Consider the residues to a given 
modulus of numbers of the form x*, where k takes the successive 
ΔΝ a ae 

The powers of x increase fairly rapidly as k increases but as 
noted before, in finite arithmetic it is not necessary to calculate these 
out in full to obtain the required residues. Thus the residues (mod 5) 
to 2*, thatis to 1, 2, 4, 8, 16, 32, 64,. . ., are obtained simply by 
multiplying each fresh residue as it is found by 2, discarding of 
course all multiples of the modulus from the product. In the present 
case we thus arrive at: 


1,2x1=2,2x2=4,2x4=8 and 
8—5=3,2x3-5=1 
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The residues (underlined) are then 1, 2, 4, 3, and the cycle starts 
again. In other words the Ath power residues of 2 to the modulus 
5 are 1, 2, 3, 4. 5 is of course a prime but do not be hasty in 
assuming from this that the residues modulo p always contain the 
integers 1,2, 3,. . .(p — 1); those of 2* (mod 7) for instance are 
1,2, and 4, only. 

The following tables (16a, δ, c) will serve to illustrate some of the 
properties of the Ath power residues; they list the residues to the Table 16b 
prime modulo > 4 < 24 of the integers 2, 3, and 10 with an 


extension in the case of 16c. Residues to 3} 
Table 16a kK QGned)s 7 tw it BB 
0 Ar ee RS a oe 
Residues to 2” 1 ιν». Gt «νιν a 
2 δ Ὁ. 2 2 2 Ὁ 
K (ὐσῶ 5 7.11 1 45 25 3 26.3 43.8 Β @ 
0 ΕΓ a es a ee 4 δ ἡ is. .3 ΝΣ 
= 5 1 35 88 
1 ΣΙ Δ, Ee ee Cae . ; ieee 
2 cae 2 Oo ee 
7 | ae aoe: 
3 oe ee ee ee. Ὰ ne 
4 I 5. ww we Ὁ ; a an ae 
5 2. ἃ 19 ΔΉ, i388 10 ee’ 
6 25. δ 9 Be ee 1 ig hag 
7 7 Δ, 9. 446 58 
12 4 lil 
8 ae oo Fee "Ὁ 
13 12 14 
9 6,5 18 6 14 ἘΣ 
10 1 10 av. 22 
15 6 12 
11 7 δ᾽ | 
16 Ἐ 77 
12 1 11 
17 13 
13 3 Ν : 
14 6 
15 12 
16 3 
17 10 
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Table 16c 
Residues to 10* 

k (Mod)7 11 13 17 19 23 31 37 41 43 
0 ee eee εν 1... ἡ 
1 3 10 10 10 10 10 10 10 10 10 
2 ἃ: ἢ Re Ss 7 BU ΝΗ 
3 6 ΣῊ & 1 ied 
4 4 2 ἃ BH. ἤὲ 37 24 
5 5 4. ὁ 3 Ὁ 1 25 
6 1 " © 24 6 ἢ 35 
7 5 15 14 20 6 
8 ἰδ 7 ἃ # 17 
9 7 18 wu 41 
10 4" one 3 23 
11 3 14 22 19 15 
12 a UF. ὦ 21 
13 11 13 15 9 38 
14 Β 16 22° 28 36 
15 uw @¢ Ss 1 16 
16 i ee 31 
17 a) 4 9 
18 δ 4 
19 21 40 
20 3 13 
21 7 1 
22 Ι 


* The column (mod 29) has been excluded in order to keep the table rea- 
sonably compact (it runs to the full length of 28 integers). 
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The following observations are not intended to replace the proper 
theoretical approach to this subject—this is competently dealt 
with in a number of excellent textbooks—but their object is to 
whet the appetite and to show some examples of how different 
branches of arithmetic interlock, sometimes in the most unexpected 
ways. 

Probably the first thing that will be noticed in these tables is that 
whilst there obviously cannot be more than (# — 1) residues to any 
given modulus π, the actual number of residues is in many cases 
only a fraction of (ἡ — 1). What is more, whenever the kth power 
residues are only ἐ(η — 1) in number they are exactly the same 
values as the quadratic residues to the modulus n. (See Table 15.) 
Looking further into this it will be seen that when there is a full 
complement of (n — 1) residues to any modulus, those of the even 
indices, (k = 0, 2, 4, . . .) are quadratic residues and therefore 
those of the odd indices are non-residues. 

Since in these tables we are dealing with residues to the powers 
of integers—of the index k—this enables us to apply to the 
multiplication of residues a logarithmic procedure similar to that 
used in ordinary arithmetic. Thus, to find the product (mod n) of 
two or more residues it is only necessary to add their corresponding 
k values together and then find the residue in line with this new k. 

The procedure will be illustrated by taking an example from each 
table. First, from Table 16a, if we take two values in the (mod 13) 
column, say 8 and 12 the product of which is 96. Then 96 = 5 
(mod 13) = 2%. The k values corresponding to ὃ and 12 are 3 
and 6, which added together equal 9. 

Similarly, (Table 165) to the modulus 11: 


3x9=27=5= 3%, and 1 +2 =3. 
Again, (Table 16c) to the modulus 19: 
6 x 17 = 102 =7 = 1013, and 4+ 8 = 12. 


(6) It is now possible to provide an explanation for the rule given 
previously concerning the multiplication of quadratic and non- 
quadratic residues. The addition of two even or two odd k’s produces 
aneven k, and hence a quadratic residue, whilst the addition of one 
even and one odd & is always an odd value and therefore a quadratic 
non-residue, 
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Attentive readers will by now, I am sure, have noticed that the 
sequence of residues (mod 19) to 2* and 10* are identical but in 
reverse; how many, I wonder, have associated this sequence with 
that described in earlier chapters (a) for determining whether a given 
number is divisible by 19, and (δ) as part of one process for con- 
structing the recurring decimal of 1/19. 


EXERCISES 


1. Show that 101 is a divisor of 99999999, 

2. Show that 252 — 1 is divisible by 23. 

3. The following five numbers all have the correct ‘endings’ and 
digital roots of square numbers. Given that the quadratic 


residues of 7 and 11 are respectively (0, 1, 2, 4) and (0, 
1, 3, 4, 5, 9) find two of the five which cannot be square. 


a 2544025 
b 2908456 
c 2712609 
d 2893401 
e 2354464 


2 FACTORISATION 


(1) A large part of that branch of mathematics somewhat loosely 
termed ‘Number Theory’ is either devoted to, or has its origins 
in the quest for prime numbers. The determination of whether a 
given (large) number is prime, or if not what are its prime factors, 
is frequently a matter of interest and more often than not one of 
considerable difficulty. In this chapter we shall examine some 
of the methods which have been developed for reducing these diffi- 
culties. 

Large numbers or small, there are of course some elementary 
observations which can be made at once. Even numbers and multiples 
of 5 are recognisable at sight, multiples of 3 respond to a simple 
and rapid mental check and thanks to a lucky combination of 
factors those of 7, 11, 13 and 37 are found with little effort. 

As we have seen in an earlier chapter many other of the smaller 
primes can be eliminated without having to use the laborious process 
of testing by actual division. These methods obviously dispose of 
vast quantities of the natural numbers—even numbers and multiples 
of 3 alone take care of two-thirds of them, for instance—but there 
remains an infinity of integers which are either prime or have 
larger prime factors than the methods used in chapters 1 and 2 are 
equipped to deal with. These can be divided broadly into two classes; 
in the one case numbers of a particular structure such as, say, those 
expressed generally by the formula x" +1 for which the ‘form’ 
of possible factors can usually be determined fairly easily. On the 
other hand there are amorphous numbers, or shall we say numbers of 
no known histury, whose factor forms can only be found by con- 
gruence techniques which for numbers of six or seven digits are 
inferior to other methods of factorisation and are quite impracticable 
for larger numbers.* Incidentally, it must be noted that in testing a 
given number WN for possible prime factors it is not necessary to 

* See Advanced Algebra, Barnard & Child, pp. 192-3. 
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try a divisor greater than +/N_ since if there is a factor larger than 
this there must also be a smaller one which has already been disclosed. 

Since testing by direct division by the primes in sequence is 
apparently only really applicable to numbers already well covered 
by published factor tables it would seem that there is little point 
in pursuing this method further. As we shall see later, however, 
in one of the most successful methods of factorisation it is of the 
utmost importance to determine the highest possible limit of primes 
which are not factors of the number we are trying to crack. Apart 
from this it will be obvious that if we have eliminated all the primes 
less than the cube root of a given number WN then there can be at 
most two factors of N and one of these must lie between ὃ. Ν 
and 4/N, clearly a piece of valuable information. 


(2) The best way of testing a number of small prime divisors ‘en 
bloc’ is to employ Euclid’s Algorithm for finding the H.C.F. (highest 
common factor) of two numbers. In essence this depends upon the 
fact that if two numbers have a factor in common then the remainder 
after dividing one by the other will also contain this factor. Similarly, 
on dividing this remainder into the previous divisor if there is a 
further remainder this will also contain the factor, and so on. 
Continuing in this manner the position is finally arrived at where 
either the remainder ry, equals 1, in which case both numbers are 
relatively prime, or rz-1 is a multiple of r,. In the latter event the 
integer ra is the greatest common factor of the two numbers. 

As so often happens in arithmetic the method is far simpler to 
operate than to describe and an elementary example should clarify 
the above explanation. Taking the numbers 21 and 56, we proceed 
to divide thus: 


56/21 = 2, + 14; 21/14 = 1, + 7; 14/7 = 2 exactly, 


and therefore the last divisor—7—is the H.C.F. of the two numbers 
21 and 56. 

It will be obvious then from the above that if we prepare a com- 
posite number containing as factors all the prime divisors we wish 
to test as possible factors of a given number N, Euclid’s Algorithm 
will supply the answer. It is conventional to express the product of 
successive primes by the Greek capital II, and specifically for 
instance the product of 3. 7. 11. 13... . . 101 is written II(p). (For 

8<p<l0l 


LL ΕΞ ΞΣ.... : I II sl . -.--- «- = 
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obvious reasons 5 need not usually be included.) The following 
example will show the basic principles of this method of factorisa- 
tion. 


Let N = 943. 


Now 943 is clearly not a multiple of 3, nor yet of 7, 11, or 13. 
And since it is less than 31? its only possible factors are then 
17, 19, 23, or 29. The product of these four primes is 215441. 
The algorithm then proceeds; 


943)215441 
1886 
~ 2684 
1886 
“7981 
7544 
Remainder = "437/943 
874 
Remainder = ~ 69)437 
414 
Remainder = ~-23)69 
69 
“0 


And therefore 23 is a factor of 943. 

The following illustration is given, not with the intention of boring 
the reader, who can easily skip it if he feels inclined, but in order to 
give those interested some feel of the ‘weight’ of the calculation for 
numbers of a moderate size. The ‘time’ advantage over trial divisions 
by the separate primes increases with the size of N. 

We shall examine the number N = 20,699,411, and since the 
primes 3, 7, 11, 13, and 37 can be so readily eliminated it will be 
assumed that they have been already tested and we then formulate 


II(p) (ex 37) = 29330271959743. 
3 


17<p<5 


The algorithm then proceeds: 


2069941 1/29330271959743 


20699411 
86308609 
82797644 
35109655 
20699411 
144102449 
124196466 
199059837 
186294699 
127651384 
124196466 
34549183 
20699411 
13849772/2069941 1 
13849772 
6849639/13849772 
13699278 
150494/6849639 
601976 
829879 
752470 
77409/150494 
77409 
73085 
73085/77409 
73085 
-4324/73085 
4324 
29843 
25944 
3901/4324 
3901 
~ 423/3901 
3807 
~ 94/423 
376 
47/94 


94 
0 


Factorisation 
Therefore 47 is a factor of N. 
Dividing out, we have N = 47 x 440413 


To continue with the next few possible prime divisors we take 
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59.61.67.....83 = Π(ρΡῚ = 8472681192817. 
59<p<83 
440413/8472681192817 
440413 
4068551 
3963717 
1048341 
880826 
1675159 
1321239 
3539202 
3523304 
1589881 
1321239 
2686427 
2642478 
43949/440413 
43949 
923/43949 
3692 
7029 
6461 
568/923 
568 
355/568 
355 
213/355 
213 
142/213 
142 
71/142 
142 


0 
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Therefore 71 is a factor of 440413. 
Thus N = 20699411 = 47 x 71 x 6203. Moreover since ν 6203 is 
less than 80 and all primes up to and including 83 have been 
tested then 6203 is prime and we have the complete factorisation. 

For dealing with large numbers the time consumed by this proce- 
dure is relatively trivial and it is generally preferable to employ 
stages of much larger II(p)s. 

To save some duplication of effort I give below the products 
of primes which I have found most convenient in practice. 


Table 17 
Products of Successive Primes II(p) 


Πρ) 

17 <p < 89 (ex 37) 
97 <p <173 
179 <p <251 
257 < p < 337 
347 <p < 419 
421 <p < 499 


21391889098600094293679777521 
700878088 17907376609 18902153322941 
385775665012409 154140882380977413 
231050535484280435538943 196422625393 
3236374326887641816997742721494821 
17861723378516715340736738647007161271 


There is no need to be intimidated by the size of these numbers; 
the first division may take a little time but after the first remainder 
is found the rest of the algorithm requires only a matter of minutes. 
Generally speaking a number of about ten or twelve digits can be 
tested for all the prime divisors less than 500 in under the hour and 
if a hand or desk calculator is available this time will be considerably 
reduced. 

Naturally if we are dealing with a number whose only possible 
factors are of a particular ‘form’ it is only necessary to use primes of 
this form in preparing multiple products for use in the algorithm. 
This means that either the time involved can be cut or the tests can 
be carried to much higher values of p. 


(3) We turn now to a method of factorisation which is infallible, at 
least in theory. The qualification must be made because there are, 
of course, random numbers of such magnitude that the labour 
involved in any computational method is prohibitive. 

This method is practicable for numbers of up to about twenty 
digits, although it must be admitted that if the number turns out to 
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be prime, as in the case of 719, the amount of man-hours required 
for the test calls for some dedication. 

The procedure is based on an identity known to every schoolboy, 
namely that 


x2 — yp? = (x + yx — y), x> py. 


It is easily shown that all odd numbers can be expressed as the dif- 
ference of two squares. For if N, an odd number, is equal to ab, 
then a and ὁ are both odd and (a + δὴ) and (a — δ) are there- 
fore even. Then 


(“$*)'- (*5") -F-w=w 


2 2 


When NV is prime then x + y=WN, and x—y= 1, as seen 
in the example 13 = (7 + 6)(7 — 6) = 72 — 65. Composite numbers 
can also be expressed in this manner but in addition, in at least one 
other way. Thus when N = 105 we have: 


105 = 532 — 522 = (53 + 52)(53 — 52) = 105 x 1 
= 192 — 162 = (19 + 16)(19 — 16) = 35 x3 
= 132— δὲ (13... 8\(13— 8)= 21x 5 
=112— 4@=(11+ 4(11— 4)= 15Χ7. 


In order then to divide a given odd number N into two factors it 
is necessary to find two squares which have WN as their difference. 
In practice this is effected by subtracting N in turn from a rising 
sequence of square numbers until a square turns up in the remainders. 
To take a simple example, let N = 217. The first square greater 
than 217 is 153 = 225, and we might then proceed as follows: 


153 -ο 225 162=256 172=289 182=324 192 = 361 
ΝΣ ΤῊ πὴ ὴΡ 
8 39 72 107 114 


Now 144 = 122, and therefore 217 = 192 — 122 = 31 x 7. 

This example serves to illustrate the basic principle of the method 
but the technique is far too clumsy to apply to large numbers. How- 
ever, a simple manoeuvre enables us to set out the work in a much 
more compact and workmanlike fashion. 

Since (n + 1)? = n? + (2n + 1), successive squares increase by 
successive odd number increments and in consequence the dif- 
ferences between the rising squares and N will increase by the 
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same increments. Thus the above individual subtractions can be 
replaced by cumulative sums in the following manner. 


152 = 225 
N =217 
ae 
(2x 15)+1= 31 
162—N= 39 
33 
nat ἜΣ 
35 
etc. 107 
37 
144 = 122 


(4) The above introduction to this subject will be expanded later, 
but first it is obviously important to be able to recognise a square 
when one is found. When dealing with moderately small numbers it is 
only necessary to have at hand a table of square numbers of which 
many are available. As N becomes large, however, one must 
eventually resort to the traditional method of square root extraction. 
Even so there are some observations to be made which will reduce 
this necessity to a minimum. 

It has already been noted that the digital roots of square num- 
bers—that is, the quadratic residues of 9—can only be 1,4,7, and 9 
and this is a test which can be applied with speed. Examination of a 
table of squares will also show that the last two or three digits of 
square numbers are confined to certain definite patterns. Using the 
convention that 0 represents any even digit and 1, any odd digit, 
it will be seen that the only possible ‘square endings’ are 01, 04, 16, 
09, 00, and 125, 225, 625. A still closer inspection reveals that a 
further restriction can put upon numbers with the endings 01 and 
09, namely that squares in these categories can only have the follow- 
ing forms of endings: 


001 - 009 
121 129 
041 . . 049 
161 169 


== 


i ee 


5 ee. ᾿ S/F: —a 
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Numbers satisfying these conditions are still not necessarily 
square but if the suspect number is large it may be preferable to make 
further tests by division by small primes rather than resort to the 
tedious process of root extraction. For instance if on division by 
7 the remainder (residue) is not 0, 1, 2, or 4 (see Table 15), the 
number is not square. 

Taken together the above tests will eliminate a considerable num- 
ber of integers which cannot possibly be squares. 


(5) Weare now equipped to attempt the factorisation of a moderately 
large number and I have chosen one of sufficient magnitude to 
demonstrate both the power of the method and also to leave some 
elbow room for the development of some useful short-circuiting 
techniques. 

Let, for example, N = 1,002,387,143. To begin with the smallest 
square greater than N is found by the customary square root 
extraction process, thus: 


3166 1? 
1002387143 
9 
61 102 
61 
626 4138 
3756 
6326 38271 
37956 


63321 31543 (to ensure that the 
63321 square is larger than 
—31778 N) 


We now have N = 316612 — 31778 and the successive differences 
from the rising sequence of squares are found as follows: 
(Note that 2 x 31661 + 1 = 63323) 


AS A Te 6 ὦὕὖὦ 
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31778 
63323 
95101* = 316622 — N 
63325 


158426 
63327 


221753 
63329 
285082 
63331 


348413 
63333 


411746 
63335 

475081* = 316682 — N 
63337 


538418 
63339 


601757 
63341 
665098 
63343 


728441* = 316722 — N 
63345 


791786 
etc. 


It will be seen (the reader can confirm this by continuing the 
column) that after ten steps the final digits of the sums begin to 
repeat themselves in the same order, and that in addition the penulti- 
mate digits have the same parity as their predecessors. Furthermore 
we note that in the intervals of ten steps only two sums have the two 
terminal digits of a possible square*. These correspond with the 
numbers 316622— N, and 31668?— N and are 95101 and 
475081 respectively. 

Owing to the cyclic pattern of the ‘endings’ it follows that the only 
sums in which squares can possibly appear are those related to the 
squares 316627, 316722, 316827,... and 316687, 316787,. .. 


inant 


ES “--“-ἰἴ - a 
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etc. We can therefore dispense with the intermediate steps and set 
up two columns, each of which will advance at ten times the rate we 
started with. 

This presents no difficulties since we are using here an Arithmetical 
Progression and the total sum of ten such steps is readily found by 
coupling them in five equal pairs, thus: 


63325 63327 63329 
63343 63341 63339 


126668 126668 126668 etc. and 5 x 126668 = 633340. 


(Note. Each successive ‘ten steps’ increment increases by 200, 

and in fact it is easily shown by A.P. formulae that when, as in this 

case, the common difference per ‘step’ is 2, then the incremental 

increase for n steps is 2n?. Thus the increase for one step is 2, 

for 10 steps—200, 20 steps—800, 90 steps—16200, and so on.) 
To proceed we now set up the following columns: 


95101 = 316622 — N * 475081 = 316682 — N 


633340 633460 
* 728441 = 316722 — N 1108541 = 31678? — N 
633540 633660 
1361981 ete. *1742201 etc. 
633740 633860 
*1995721 2376061 
633940 634060 
2629661 *3010121 


None of these sums are squares but we now note that whilst they 
all have the required last two digits of squares only those starred have 
the necessary parity for the third digit from the last. Successive A.P. 
numbers are always rigid pattern-followers and if the above columns 
are extended it will still be found that only every other sum has the 
three permissible terminal digits of a possible square. 

It is therefore possible to reduce the number of successive additions 
still further, this time to cover twenty of the original steps. As the 
following columns will now have all the correct terminal digits of 
squares we shall pay attention to the digital roots in the next summa- 
tions. Starting with the starred members of the last columns we then 
proceed: 


118 Factorisation 
D.R. D.R. 
728441 = 316722 — N 8 475081 = 316682 — N 7 
1267280 1267120 
1995721 = 316922 — N 7 1742201 8 
1268080 1267920 
3263801 5 3010121 8 
1268880 1268720 
4532681 2 4278841 = 317282 — N 7 
1269680 1269520 
5802361 = 317527 — N 7 5548361 5 
1270320 
6818681 2 
1271120 


8089801 = 317882 — N 7 


Note that as each of the above additions is equivalent to twenty 
of the original steps the increments increase each time by 800. 
It will also be seen that the only digital root of a square to appear 
in this cycle—namely 7—occurs at every third sum (i.e. at every 
sixty of the original steps). Consequently there is no point in looking 
for squares in the sums of anything but the (31692 + 60k)? — Ν, 
and (31668 + 60k)? -- N cycles, and as no squares have yet 
appeared we now set up two further columns which in addition to 
having the required terminal digits will contain only sums possessing 
the digital root of 7. (Incidentally, since each of these sums will 
now correspond to 60 of the original steps each addition increases 
by 2 x 60? = 7200.) 


1995721 = 316922 — N 4278841 = 317282 — N 
3806640 3810960 


5802361 = 317522 — N 8089801 = 317882 — N 
3813840 


9616201 = 318122 — N 
We now have a square at last, for 9616201 = 31012. 
Hence N = 1002387143 = 31812? — 31012 


= (31812 + 3101)(31812 — 3101) 
= 34913 x 28711. Which is the complete factorisation, 
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Shorn of explanatory verbiage and repetitions in the columns the 
above process can be seen to be remarkably economical in both time 
and paper: it is still open to further refinement as will be shown later. 
At the same time modifications of the technique may sometimes be 
required; with some Ns, for instance, the (square) digital roots 
are not evenly spaced and this entails setting up four columns at 
this stage. There should be no difficulty, however, in extemporising 
to meet any situation if the above arguments have been thoroughly 
mastered. 


(6) A fundamental difficulty arises in this method of factorisation 
when WN has only two factors which have widely different orders of 
magnitude. When this happens, or in the extreme case when JN is 
both large and prime, the operation becomes prolonged, perhaps to 
an intolerable extent. 

If the process then is to remain practicable in such circumstances 
it is essential to find some way of at least reducing this difficulty. 

Now it will be clear that however many factors N may possess 
the method described above produces only two; one or both of 
these may of course be composite. 

Let N = ab, and as we want N in the form x? — y? then we 
must have a= x+y, and ὃ τὰ -- γ. 

(It should be noted that since N is odd—or can be made so by 
division by some power of 2—then both x + y and x — y must 
also be odd.) 

If a and ὃ are of the same order of magnitude we can say that 
a~b~-+/N, which implies that y is relatively small and will be 
found early in the proceedings. On the other hand if a and ὁ differ 
appreciably then it can be said that a= (approximately) kb, 
where &k is an integer. Obviously if k can be found then kN will 
have the two nearly equal factors a and kb which will quickly 
be disclosed. 

Unfortunately there is no way of predicting the value of k and 
so it can only be found by trial. That is, if the above algorithm 
shows no signs of finding a square it should be abandoned and a new 
test started on 3N, then on 5N, and so on. (The reason why even 
values of k are not recommended can be explained as follows. 
When & is odd then both x + y and x — y must be odd and the 
condition noted above is satisfied. When, however a is nearly 
equal say to 2b and the algorithm is applied to 2N no solution will 
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emerge because we have now made one of the factors even whilst 
the other remains odd and it is no longer possible for them to assume 
the forms x + y and x — y. Furthermore, if 4N is tried, although 
now both new factors will be even they will be in the same ratio as 
before and no advantage will have been gained. In such a case, 
therefore, N must be multiplied by 8 before the new factors 
become nearly equal, and in general if k is an even number N 
must be multiplied by 4k.) 

To recapitulate, it has been shown that if N = ab and a ~kb, 
then the difference of squares method will quickly produce the 
factors of N if it is applied to 


kN when k is an odd number, 
4kN when k is an even number. 


It may come as a surprise that increasing the size of the number to 
be tested actually leads to a quicker solution of the problem but 
the reader can easily satisfy himself of the truth of this. 

It will now be seen why it is important to eliminate as many of the 
smaller prime divisors as possible before beginning this procedure 
since we are then able to define the lowest limit of k. Thus when it is 
clear that there are no factors less than, say W, the first k to be 
tried should be just greater than N/W?. 


(7) The foregoing process is powerful and yet requires little more than 
an attention to detail which one must accept in this brand of mathe- 
matics. After a few practice shots at numbers of say, three to six 
digitt—Premium Bond numbers provide excellent material—the 
reader should be able to tackle quite large numbers with confidence. 
In fact he will be equipped to perform what some writers have re- 
garded as near-miracles. 

On the subject of primes, and referring to an often quoted Fermat 
legend (repeated even as recently as March 1964 in Scientific 
American in almost identical terms), Kasner and Newman in 
their stimulating book Mathematics and the Imagination have this 


to say (p. 187): 


‘Curiously enough there is reason to believe that certain mathe- 
maticians of the seventeenth century, who spent a great deal of 
time on number theory, had means of recognising primes unknown 
tous,. . .itis still a source of wonder that Fermat replied without 
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a moments hesitation to a letter which asked whether 100895598169 
was a prime, that it was the product of 898423 and 112303, and 
that each of these numbers was prime. Without a general formula 
for all primes, a mathematician, even today, might spend years 
hunting for the correct answer.’ 


The ‘difference of squares’ method, which was of course well 
known to Fermat, breaks this number down almost before we can 
get started, but look what happens when we apply the process to 
8N. I give the calculation in full in order to show that the time in- 
volved is trivial. 

We have WN = 1008955981169; 8N = 807164785352. Then 
extracting the square root of 8N: 


807164785352(898424 
64 
169 1671 
1521 
1788 ~ 15064 
14304 
17964 ~ 76078 
71856 
179682 "422253 
359364 
1796844 6288952 
7187376 
— 898424 
Thus ΒΝ = 8984242 — 898424 


= 898424 (898424 — 1) 
= 8 x 112303 x 898423 


Therefore N = 112303 x 898423. 


So much for ‘lost methods’, ‘general formulae for primes’ and 
‘years of hunting’. 


(8) Before ending this chapter it will perhaps be appropriate to 
include a note on ‘checking the work’. In calculations of this kind, 
which may on occasion become protracted, it is always well at 
intervals to check the accuracy and the progress of the operation. 


Ee 


122 Factorisation 
This is easily carried out and it will be enough to give only a simple 
example. 

Suppose the ‘difference of squares’ algorithm is being applied to 
N = 3527, a number which is easily seen to have none of the 
factors 3, 7,11, and 13. The working starts with 


N = 60? — 73 
and proceeds to the point 
N = 126? — 12349 


without a square number appearing. At this stage it is decided to 
make a check. 


The largest square less than 12349 is 1112 = 12321 which differs 
from 12349 by 28. 
Now: 


1262 — 1112 = 15 x 237 = 3555, 
And 


3555 — N = 28. 


Therefore the calculations are without error and furthermore the 
check shows that all the primes between 4/N and 15 have been 
tested. It follows that 3527 must therefore be prime. 

This method of checking should never be neglected since it not 
only confirms the accuracy of the work but at the same time provides 
an indication of when one should start again with a new value of KN. 


EXERCISES 


1. Find one factor common to both 16733 and 35699. 
2. N = 1181027 has one factor less than 30. What is it? 


3. Using the ‘difference of squares’ method find the factors of 
N = 1093709. 


4 3 MORE FACTORISATION 
METHODS 


(1) We are all familiar with the ‘square’ numbers, which can be 
represented graphically in the form: 


ἙΝ 
4 


9 


Ι6 εἰς. 


‘Triangular’ numbers are those formed by the cumulative sums of the 
rows in a triangular array, thus: 


* 1 

ok tk 3 

* ee * 6 

* * Ἕ Ἂς 10 

τς ΑΝ ΜΠ 15 etc. 
and they produce a sequence of which the th term is 

n(n + 1) 
ΣῈ 


These numbers share with the squares many useful properties, 
not all of which seem to have received full recognition. 

In my opinion their employment in factorisation techniques 
is vastly to be preferred to that of the square numbers. It is easily 
seen that, just as with the squares, any integer which can be expressed 
as the difference of two triangular numbers is composite and, when 
so resolved, the factors can be found immediately. 
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Thus, (x > y): First find the least triangular number 5» ΝΟ by extracting the 
square root of 2N. 
χα ἘΔ γυ 1 x*-—y?+x-y 
ΒΕ ἈΠ ARETE Aah at 23 ἐ 
2N = 875782 
_ =e tyro —%) 81 
2 13 6057 
549 
= ἀπ ΣῈ 1866 10882 
11196 
When using these numbers for factorisation the same use can be 936 x 937 
made of digital roots, terminal digits and multiples of N as those and we have —————. = 438516. (= say, 936) 
described in some detail in the last chapter, with, of course due 2 


recognition of the basic difference in these properties. Thus, tri- 
angular numbers have digital roots of 1, 3, 6, and 9, and can 
only end in 0, 1, 03, 53, 5, 6, 28, 78. 

In this system the ‘step’ increments increase by 1, instead of 2, 
but the digital roots, etc., have exactly similar cyclic properties and 
offer the same opportunities for abbreviation. The method has one 
disadvantage, namely that tables of triangular numbers are not so 
readily available as those of squares. They are, however, easily 
constructed particularly if one has access to a hand or desk calculat- 
ing device, and in any case the loss is not all that important when 
dealing with large numbers. A suspected triangular number, that is 
one having the right ending and digital root, is tested by first 
multiplying by two and then extracting the square root. Thus: 
T = 11935. Then 11935 x 2 = 23870. 


15 4 
2T = 23870 = 154? + 154 = 154 x 155 
Ι and therefore T = (154 x 155)/2 
5 158 which is of triangular form 
125 
304 1370 
1216 
154 


After what has been said there should now be no difficulty in 
understanding the working of the following example. 
To find the factors of N = 437891 (= 397 x 1103). 
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The calculation then continues: 


438516 
—N 437891 


625 = 936; — N 
937 


4379 (940, — N) 


11935 = (948, — N) = 154 


We now have, N = 948; — 154; and therefore N is composite, 
its factors being found as follows: 


948 948 
+ 155 — 154 
1103 794 = 2 x 397 


It will be noticed that as each value of (x; — N) is found x is 
numerically equal to the previous incremental addition; because of 
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this it is simpler to keep track of the operation than in the difference 
of squares method. The greatest advantage that the triangular num- 
bers have over the squares in this system of factorisation lies in the 
number of steps required for ultimate solution. In the above example 
the factors were found after only twelve steps whereas the “difference 
of squares’ process requires 87 steps for the same N. 


(2) We turn now to some less orthodox systems which can be used 
for finding factors. They are included more for their individual 
interest than for their practical value in attacking large numbers, 
although any of them might well repay some development work. 

The first is simply a variant on the method described in the last 
chapter. It is deliberately being treated separately because unlike 
the former method its power is seen to the best advantage when WN is 
either prime or has factors of widely different orders of magnitude— 
that is, when k proves to be large. 

The principle becomes clear when it is remembered that whilst 
there are exactly p—1 quadratic residues, all of them different, 
to any prime p, the number of residues to a composite number M 
is always less than M — 1 and some of them are repeated at least 
once. 

Briefly, the process starts with a progression of steps produced 
in exactly the same way as in the ‘difference of squares’ method, 
but in this case the progressive totals are not allowed to exceed N 
in value. When this occurs N is subtracted and the steps are then 
continued in the normal way. From this point onward a watch is 
then kept for any sum which has appeared before and as soon as a 
duplication is seen the factorisation can be completed. 

To give an example of the process, we have: 


Ν = 2533 = 512 — 68. (= 17x 149. k~9) 


128 More Factorisation Methods 
Then, 
68 1311 118 
103 125 145 
171 1436 263 
105 127 147 
276 1563 410 
107 129 149 
383 1692 (652) 559 (752) 
109 131 151 
492 (552) 1823 (662)* 710 
111 133 153 
603 1956 863 
113 135 155 
716 2091 1018 
115 137 157 
831 2228 1175 
117 139 159 
948 2367 1334 (802) 
119 141 161 
1067 (602) 2508 1495 
121 143 163 
1188 2651 > Ν 1658 
123 2533 165 
1311 18 1823 (832)* 
* The sum 1823 has been repeated and the factors are found thus: 
N = 662 — 1823 
2N = 832 — 1823 
Subtracting N = 832 — 662 = (83 + 66)(83 — 66) 
And then N = 149 x 17. 


As a suggestion for further research into this subject, it might be of 
interest to note that the residues to triangular numbers have similar 
properties to quadratic residues and are also duplicated, sometimes 
at frequent intervals, when N is composite. 


(3) The following algorithm is not so much a method of factorisa- 
tion as a device for reducing the size of the number under test. 
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I have not seen it mentioned before but that, of course, does not give 
it any claim to originality. 

To find the factors of N first set up a column (A) in which the 
first term is (N —1)/2, the following lines being successively 
reduced by unity. At the same time odd integers 1, 3, 5, 7,... 
are set alongside these terms in a column (B). Then those numbers 
in col. B which divide their companions in col. A are also divisors 
of N. 

To illustrate, we take N = 1001. Then (N — 1)/2 = 500. 

Then proceed, 


A B 
500 1 
499 3 
498 5 
497 7 


And since 497/7 = 71, then 7 isa factor of 1001. Furthermore 
as (2 x 71) +1= 143 = 11 x 13, we have the other factors. 
Or alternatively the columns could be continued thus: 


496 9 
495 11 (495 = 11 x 45) 
494 13 (494 = 13 x 38) 


Since the integers in col. A decrease by exactly half the increase 
of those in col. B it is a simple matter to arrange for composite num- 
bers to be omitted from col. B, or to skip as many numbers as we 
like. In the first instance column B can, and indeed should, contain 
only the odd primes, and in the second whenever the possible factors 
of N are known to be of a particular ‘form’ only primes of this 
form need be tabulated. 

Thus, for example if N = 11111, and it is known, as we saw 
earlier, that only numbers of the form 305 - 1 and 30” + 1] 
are eligible then we can proceed as follows: 


(N — 1)/2 = 5555. 


Then we have 5555 1 
5550 11 
5540 31 
5535 41 
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On testing it is found that 5535/41 = 135. 

And since (135 x 2) + 1 = 271, we have the factors 41 and 271. 

This method has the following advantages: it allows little oppor- 
tunity for errors in calculation; with the aid of factor tables it is 
effective in quickly eliminating quantities of small primes particularly 
when they are of known ‘form’; and as it starts with N/2 it therefore 
virtually doubles the range of the available factor tables. 


(4) The method of factorisation about to be described is taken 
from L. E. Dickson’s History of the Theory of Numbers Vol. 1 
and therein is attributed to D. Biddle (1911). 

It is included here not for its practicability but to provide a basis for 
speculation on the devious processes which led up to its conception. 

Process. Express N as 55 - A where S® is the largest square 
<N. Write three rows of numbers, the first beginning with A, 
(or A—S if A> 57); the second beginning with S (or S +1, 
in the latter case), and increasing by 1, the third beginning with S 
and decreasing by 1. 

Let An, Bn, Cn be the mth elements in the respective rows. 
Then Cy =Ca-1—-1; Ba = Ba-1+1; An = An-1 + Ba-i — 
KCn. When (An-1 + Bn-1) is greater than Cp, a multiple of Ομ 
is subtracted so as to leave a positive remainder and then By, = 
Bn-1 ἘΚ. When a value of n is reached for which A, = 0, 
then N = Bn x Cn. 

For an example we take N = 589 = 245 + 13 


A 13 14 17 1 9 0 
B ΠῚ ἈΠ DS 
Cc ΝΣ δ᾽ Δ ΠΣ ee ἂν 


The factors of N = 589 are thus disclosed as 19 and 31. 

I am sorry if the above account lacks a certain amount of elegance 
but am comforted by Dickson’s laconic note ‘It may be best to start 
with 2N.’ 
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EXERCISES 


1, Find the divisors of the following numbers using any method of 
factorisation. 


a 15871 
b 15853 
c 15863 
2. The factors of N = 13333 are of the form 22k + 1. What are 


they ? 
3. The factors of N = 548497 are of the form 26k + 1. What are 
they? 


i A NOTE ON THE CONVERSE 
OF FERMAT’S THEOREM 


(1) Many readers must have been puzzled on reading that such and 
such a (large) number is composite but its factors are not known. 
One would think that the authors of books in which such statements 
appear—and there are quite a few of them in the ‘popular’ mathe- 
matics field—might have spared a few lines to explain the basic 
principle on which an observation of this sort depends. 

Let us consider Fermat’s theorem again. This states in effect, 
that if p is a prime and a is neither equal to, nor a multiple of p, 
then (a?-1— 1) is a multiple of p. Or in congruence language 
a?-1 = 1 (mod p). All this is perfectly true but one must be careful 
not to infer that the converse is necessarily true, namely that if 
aN-1 =1 (mod N) then N must be prime. The probability that 
it is, is very high but in spite of this it can be proved that there are an 
infinity of composite N’s which satisfy the congruence. 

Lucas first laid down, and in 1891 proved the necessary conditions 
for a true converse of the theorem, namely that ‘If αὖ = 1 (mod N) 
for x = N—1, butzotfor x a proper divisor of N — 1, then N 
is a prime.’ 

Since then Lehmer in particular has developed many refinements 
to the theory, so much so that it can be considered now as one of the 
few valid tests for the primality of large numbers. The arguments are 
somewhat advanced and quite beyond the scope of this book; 
probably the best presentation is to be found in D. H. Lehmer’s 
paper, “Tests for primality by the converse of Fermat’s theorem’, 
Bull. Amer. Math. Soc., vol. 33, 1927, pp. 327-340. 

On the other hand the converse presents no difficulties when applied 
as a test for composite numbers, and it is legitimate to conclude that 
if αν πὶ =r (mod N), where r proves to have any other value than 
1, then Ν᾽ is not prime and is therefore a product of two or more 
primes. 

To carry out this test is is necessary to find the remainder on 


A Note on the Converse of Fermat’s Theorem 133 


dividing aN-1 by Ν᾽ and as a can take any value prime to N 
it is convenient to let a = 2. The actual calculation uses the follow- 
ing properties of a congruence, namely that if 2 ΞΞ Καὶ (mod N) 
then 2? = R?, 2° = 2R?, etc. It is advisable to set out the operations 
in an orderly manner and the following explanation and the three 
worked examples—although applied to only small N’s—will indicate 
the procedure to be followed in testing numbers of any magnitude. 
The calculation employs two columns of figures labelled here A 
and 24 (mod WN) respectively. The left hand column (A) is con- 
structed first and starts with N — 1. Each successive line is obtained 
by dividing the previous one by 2, (after subtracting 1 if the 
number is odd). The last entry is, of course 1 and against this the 
next column is started from the bottom with the number 2. (It is as 
well at this stage to make a mark against the odd numbers in Col. A 
as a reminder in what follows.) 

The rising sequence of this second column is now formed by 
successive squaring, the resultant squares also being doubled where 
the corresponding entries in Col. A are odd. Each value is then 
divided by N and the remainder entered in the next line above. 
The final entry at the top of the column is then the residue (mod NV) 
of 2%-1! and if this is any number other than 1 then N is compo- 
site. 


Example 1 
N = 437 
A 24 (mod N) 
436 359 
218 213 
109* 173 
54 58 etc. 
27* 170 2 Χ 1112 = 24642 = 170 (mod 437) 
13* --111 2x 642 -- 8192 = 326 (mod 437) see 
below 
6 64 82 
3* 8 2x Ὁ» 
1" 2 


Therefore 2437-1 = 359 (mod 437) and N is composite. It 
should be noted that it is sometimes more convenient to use a 
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negative value as the remainder since this still becomes positive 
again after squaring. Thus 2 x 64% = 8192 = 326 (mod 437) and 
326 — 437 = -- 11. It is simpler to square —111 than 326, a 
saving of time which can be significant when the modulus is large. 


Example 2 
N = 347 

A 24 (mod N) 
346 1 
173* —1 
86 120 
43* 193 
21* 231 
10 —17 
5* 32 
2 4 
1 2 


Therefore 347 may be prime. (In fact it is.) 


Example 3 
N = 561 

A 24 (mod N) 
560 1 
280 1 
140 67 
70 166 
35* 263 
17* —202 
8 256 
4 16 
2 4 
1 2 


Although it is seen that 2559 = 1 (mod 561), 561 is not prime. 
It is in fact equal to 3 x 11 x 17. 

There is, of course, no need to confine the value of a to 2, 
though for obvious reasons it is best to keep to small values. If, 
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for any reason one chooses to use a = 3, then the base of the 
right hand column is started with 3, otherwise the procedure follows 
exactly the same pattern as with a = 2. 

Again, there is no need, as can be seen in the above examples, to 
continue the col. A downwards beyond the point where a4 becomes 


less than N provided suitable power tables are available. (See 
Appendix.) 


I 5 CONCLUSION 


(1) In what has gone before I have introduced to the reader a very 
small section of that branch of mathematics known as the Theory 
of Numbers. The main theme of the book has been factorisation, a 
subject which, once the bug has bitten is intensely fascinating and 
yet at the same time embedded in difficulties. I have hoped to 
give some indication of how some of these difficulties may be 
chipped away and to have stimulated an interest in research into 
the subject. 

Let no one be deterred from this by the advent of modern electronic 
computers and their spectacular success in finding ever larger primes 
of a particular ‘form’; the computer may take only minutes to 
provide the answer to a given problem but its actual programming 
may well have required hundreds or even thousands of skilled man- 
hours to prepare. Nor should we forget that most of these programmes 
depend for their success upon principles laboriously hammered out 
and tested by mathematical enthusiasts throughout the ages. 

A chapter on some of the known properties of primes was almost 
essential in this context. They have worried mathematicians, parti- 
cularly amateurs who presumably have more time to spend on 
matters of doubtful outcome, since their ‘separateness’ was first 
recognised. It may be that the modification to Euler’s criterion sug- 
gested in this chapter (perhaps with further necessary restrictions), 
could lead somewhere, in which case it ought not to be unduly 
difficult to extend the test to numbers of the other forms. On the 
other hand we may be asking altogether the wrong kind of questions 
about primes. Perhaps an intensive study of the more general ‘relative 
primes’ might in the long run be more profitable. Here indeed there 
could be scope for a new Einstein. 


(2) Of course mathematicians, even Number Theory addicts, do not 
confine their attentions all the time to such fundamental problems 
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and it might be fitting to close with a few samples of the many 
curious observations that have been made from time to time. 


1. Aproperty of 24; This is the largest value of n that is divisible 
by all the integers less than 4/n. 


2. A property of 30; When n = 2, 3, 4, 6, 8, 12, 18, 24, 30, 
the integers less than, and relatively prime to 7 are all unity 
and primes. No number greater than 30 has this property. 
(Uspensky & Haslet.) 


3. In the following subtractions all the original digits are repeated 
in the answers once only. 


9876543210 — 0123456789 = 9753086421 
987654321 — 123456789 = 864197532 
98754210 — 01245789 = 97508421 


4. 54+ 72+ 90 = 6° 
548 + 723 -ἰ 908 = 1088 = (72 + 90 — 54)8 


5. All the integers except the powers of 2 can be expressed as the 
sum of consecutive numbers. 


6. The addition of unity to the product of any four consecutive 
numbers produces a square number. 


7. In the number 312132 we have, 1 digit between the ls, 
2 digits between the 2s, 
3 digits between the 3s. 
Can this system be extended? 


8. Some examples of multiplications in which all the digits are 
used once only: 


7 Χ 9403 = 65821: 3 x 1458 =6 x 729. 
9. ‘Fermat’s Quotient’, (27-1 — 1)/p is only a square when 
p=3 or 7. 
10. 512 = (5 + 1 + 2)3: 
47045881000000 = (47 + 4 + 58 + 81)®. 
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11, 


12. 


13. 


14, 


15. 
16. 


17. 


18. 


19. Can the following sequence be continued, 1, 3, 8, 120,. . .? 


Conclusion 


33 +4. 48 + 53 = 63 
15 + 33 + 48 + 53 + 88 = 98 
33 4 48 + 53 + 83 + 108 = 123 
13 4+ 53 + 63 + 73 + 83 + 108 = 133 
23 + 33 + 53 + 78 + 83 + 98 + 108 = 148 
118... 12% + 133 + 143 = 208 
68 + 72 + 88+... + 688 + 698 = 1808 
11348 + 1135% + 11363 +. . . + 21328 + 21338 = 168308 
45 4+ 55 + 65 + 75 + 95 + 115 = 125 
55 + 105 + 115 + 165 + 195 + 295 = 305 
4! = 24 
5! = 120 
7! = 5040 All become squares on adding 1. 
No similar cases are known below 1020!. 


6! x 7! = 10! 


The two numbers 57321 and 60984 together contain the ten 
digits. Each of the squares of these numbers; 3285697041 
and 3719048256 contain all ten digits. 

(a) 1+2=3 

44+54+6=7+4+8 
9+10+ 114+12=13+14+15 etc. 
The highest L.H. number = πίη + 1) 
(δ) 32 + 42 = 53 
10? + 112 + 122 = 132 + 142 
212 + 222 + 232 + 242 = 252 + 265 + 27? etc. 
The highest L.H. number = (2n(n + 1))? 
But we cannot continue on these lines for 
53 + 68 = 78 — 2. 

Factors of 10" containing no 0’s. Are any to be found between 

102 = 4 x 25 


108 = 8 x 125 and 
1088 = 8589934592 x 116415321826934814453125? 


(Where the product of any two integers is x? — 1.) 


APPENDIX 


Square Roots 


(To seventeen significant figures) 


νη 
1-4142135623730955 
1-7320508075688773 
2:2360679774981797 
2-4494897427831781 
2-6457513110645906 
2:8284271247461909 
3-1622776601683793 
3-3166247903553998 
3-6055512754639893 
3-7416573867739414 
3-8729833460437668 
4-1231056256176484 
4-3657754179526917 
4-5825756931198854 
4-6904157598233230 
4-7958315233127196 
5:0990195135927848 
5°3851648071345040 


n n! 
0 0 1 3 im. 2 {δὰ 
1 1 1 S 6 δε 8. 4 te Ὧξ 
2 > 2° 20 2 253 26 200 235 88 
$ 465 496 528 S61 5305 630 66 
3 6 4 820 861 903 946 1035 1081 
: τ 5 1275 1326 1378 1431 1485 1540 [596 
6 1830 1891 1953 2016 2080 2145 2211 
6 720 7 2485 2556 2628 2701 2775 2850 2926 
7 5040 8 3240 3321 3403 3486 3570 3655 3741 
ς 40320 9 4095 4186 4278 4371 4465 4560 4656 
9 362880 10 5050 5151 5253 5356 5460 5565 5671 


13 6227020800 

14 ΒΒ7178291200 

15 1307674368000 

16 2092278988800 

17  355687428096000 

18 ὠ 6402373705728000 

19  121645100408832000 

20  2432ΘΟΖΟΟΒΙ76640000 

21 $1090942171709440000 

22 1124000727777607680000 

23 25852016738884976640000 

24 Φ620448401733239439360000 

25 [5511210043330985984000000 

26 403291461126605635584000000 
27 10888869450418352160768000000 
28 ὀ  ὠ8304888344611713860501 504ΟΟΟ000 
29 8841761993739701954543616000000 33 54615 54946 55278 55611 55945 56280 56616 
30 —-:265252859812191058636308480000000 | ἮΝ ee en ἜΝ eee Seco 
31: 8222838654177922817725562880000000 | 

33 868331761881 1886495516194401280000000 3] ens ρος GIN GTS TORS TEND TN 
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Factorial ‘n’ Triangular Numbers n(n + 1)/2 
1HRSXOMEAR. LKB n 0 1 2 3 4 5 6 
| 


34 295232799039604140847618609643520000000 39 76245 76636 77028 77421 77815 78210 78606 
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80200 
84255 
88410 
92665 
97020 


81003 
85078 
89253 


81406 
85491 
89676 
93528 93961 95266 
97903 98346 99681 


101475 101926 102378 102831 103285 103740 104196 
106030 106491 106953 107416 107880 108345 108811 
110685 111156 111628 112101 112575 113050 113526 
115440 115921 116403 116886 117370 117855 118341 
120295 120786 121278 121771 122265 122760 123256 


125250 125751 126253 126756 127260 127765 128271 
130305 130816 131328 131841 132355 132870 133386 
135460 135981 136503 137026 137550 138075 138601 
140715 141246 141778 142311 142845 143380 143916 
146070 146611 147153 147696 148240 148785 149331 


81810 
85905 
90100 
94395 
98790 


82215 
86320 
90525 
94830 
99235 


82621 
86736 
90951 


151525 152076 152628 153181 153735 154290 154846 
157080 157641 158203 158766 159330 159895 160461 
162735 163306 163878 164451 165025 165600 166176 
168490 169071 169653 170236 170820 171405 171991 
174345 174936 175528 176121 176715 177310 177906 


180300 180901 181503 182106 182710 183315 183921 
186355 186966 187578 188191 188805 189420 190036 
192510 193131 193753 194376 195000 195625 196251 
198765 199396 200028 200661 201295 201930 202566 
205120 205761 206403 207046 207690 208335 208981 


211575 212226 212878 213531 214185 214840 215496 
218130 218791 219453 220116 220780 221445 222111 
224785 225456 226128 226801 227475 228150 228826 
231540 232221 232903 233586 234270 234955 235641 
238395 239086 239778 240471 241165 241860 242556 


245350 246051 246753 247456 248160 248865 249571 
252405 253116 253828 254541 255255 255970 256686 
259560 260281 261003 261726 262450 263175 263901 
266815 267546 268278 269011 269745 270480 271216 
274170 274911 275653 276396 277140 277885 278631 


83028 83436 83845 
87153 87571 87990 
91378 91806 92235 
95703 96141 96580 
100128 100576 101025 


104653 105111 105570 
109278 109746 110215 
114003 114481 114960 
118828 119316 119805 
123753 124251 124750 


128778 129286 129795 
133903 134421 134940 
139128 139656 140185 
144453 144991 145530 
149878 150426 150975 


155403 155961 156520 
161028 161596 162165 
166753 167331 167910 
172578 173166 173755 
178503 179101 179700 


184528 185136 185745 
190653 191271 191890 
196878 197506 198135 
203203 203841 204480 
209628 210276 210925 


216153 216811 217470 
222778 223446 224115 
229503 230181 230860 
236328 237016 237705 
243253 243951 244650 


250278 250986 251695 
257403 258121 258840 
264628 265356 266085 
271953 272691 273430 
279378 280126 280875 


281625 282376 283128 283881 284635 285390 286146 286903 287661 288420 
289180 289941 290703 291466 292230 292995 293761 294528 295296 296065 
296835 297606 298378 299151 299925 300700 301476 302253 303031 303810 
304590 305371 306153 306936 307720 308505 309291 310078 310866 311655 
312445 313236 314028 314821 315615 316410 317206 318003 318801 319600 
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320400 321201 322003 322806 323610 324415 325221 326028 326836 327645 
328455 329266 330078 330891 331705 332520 333336 334153 334971 335790 
336610 337431 338253 339076 339900 340725 341551 342378 343208 344035 
344865 345696 346528 347361 348195 349030 349866 350703 351541 352380 
353220 354061 354903 355746 356590 357435 358281 359128 359976 360825 


361675 362526 363378 364231 365085 365940 366796 367653 368511 369370 
370230 371091 371953 372816 373680 374545 375411 376278 377146 378015 
378885 379756 380628 381501 382375 383250 384126 385003 385881 386760 
387640 388521 389403 390286 391170 392055 392941 393828 394716 395605 
396495 397386 398278 399171 400065 400960 401856 402753 403651 404550 


405450 406351 407253 408156 409060 409965 410871 411778 412686 413595 
414505 415416 416328 417241 418155 419070 419986 420903 421821 422740 
423660 424581 425503 426426 427350 428275 429201 430128 431056 431985 
432915 433846 434778 435711 436645 437580 438516 439453 440391 441330 
442270 443211 444153 445096 446040 446985 447931 448878 449826 450775 


451725 452676 453628 454581 455535 456490 457446 458403 459361 460320 
461280 462241 463203 464166 465130 466095 467061 468028 468998 469965 
470935 471906 472878 473851 474825 475800 476776 477753 478731 479710 
480690 481671 482653 483636 484620 485605 486591 487578 488566 489555 
490545 491536 492528 493521 494515 495510 496506 497503 498501 499500 
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A Few Specimen Primes 2" 3» 
In making up numbers for testing factorisation methods it is essential 
to employ factors known to be primes. In a book of this kind it is : 4 9 
neither necessary nor practicable to tabulate long lists of primes in ᾿ ὃ 21] 
sequence and the following selection should provide enough variety 4 16 81 
for most purposes. 5 32 243 
6 64 729 
113 401 797 2309 8803 13709 78901 107309 1000159 100004309 7 128 2187 
127 409 809 2311 8807 13711 78919 107339 1000171 100004327 8 256 6561 
131 419 811 2333 8819 13721 78929 107347 1000183 100004347 9 512 19683 
137 421 821 2339 8821 13723 78941 107351 1000187 100004363 10 1024 59049 
139 431 823 2341 8831 13729 78977 107357 1000193 100004389 11 2048 177147 
149 433 827 2347 8837 13751 78979 107377 1000199 100004393 12 4096 531441 
151 439 829 2351 8839 13757 78989 107441 1000211 100004407 13 8192 1594323 
157 443 839 2357 8849 13759 79031 107449 1000213 100004417 14 16384 4782969 
163 449 853 2371 8861 13763 79037 107453 1000231 100004449 15 32768 14348907 
167 457 857 2377 8863 13781 79039 107467 1000249 100004461 16 65536 43046721 
173 461 859 2381 8867 13789 79043 107473 1000253 100004473 17 131072 129140163 
179 463 863 2383 8887 13799 79049 107507 1000273 100004477 18 262144 387420489 
181 467 877 2389 8893 13807 79063 107509 1000289 100004501 19 524288 1162261467 
191 479 881 2393 8923 13829 79067 107533 1000291 100004503 20 1048576 3486784401 
193 487 883 2399 8929 13831 79103 107563 1000303 100004507 21 2097152 10460353203 
197 491 887 2411 8933 13841 79111 107581 1000313 100004519 22 4194304 31381059609 
199 499 907 2417 8941 13859 79133 107599 1000333 100004521 | 23 8388608 94143178827 
211 503 911 2423 8951 13873 79139 107603 1000357 100004533 24 16777216 282429536481 
223 509 919 2437 8963 13877 79147 107609 1000367 100004537 25 33554432 847288609443 
227 521 929 2441 8969 13879 79151 107617 1000381 100004549 - 26 67108864 2541865828329 
229 523 937 2447 8971 13883 79153 107621 1000393 100004551 27 134217728 7625597484987 
233 541 941 2459 8999 13901 79159 107641 1000397 100004561 28 268435456 22876792454961 
239 547 947 2467 9001 13903 79181 107647 1000403 100004563 29 536870912 68630377364883 
241 557 953 2473 9007 13907 79187 107671 1000409 100004629 30 1073741824 20589 1132094649 
251 563 967 2477 9011 13913 79193 107687 1000423 100004647 31 2147483648 617673396283947 
257 569 971 2503 9013 13921 79201 107693 1000427 100004651 32 4294967296 1853020188851841 
263 571 977 2521 9029 13931 79229 107699 1000429 100004677 33 8589934592 5559060566555523 
269 577 983 2531 9041 13933 79231 107713 1000453 100004719 34 17179869184 16677181699666569 
271 587 991 2539 9043 13951 79241 107717 1000457 100004741 35 34359738368 5003 1545098999707 
277 593 997 2543 9049 13963 79259 107719 1000507 100004813 36 68719476736 150094635296999121 
| 37: 137438953472 450283905890997363 
38 274877906944 1350851717672992089 
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549755813888 
1099511627776 
2199023255552 
4398046511104 
8796093022208 
17592186044416 
35184372088832 
70368744177664 
140737488355328 
281474976710656 
562949953421312 
1125899906842624 
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4052555153018976267 
12157665459056928801 
36472996377170786403 
109418989131512359209 
328256967394537077627 
984770902183611232881 
2954312706550833698643 
88629381 19652501095929 
265888 14358957503287787 
79766443076872509863361 
239299329230617529590083 
71789798769 1852588770249 


ANSWERS TO THE EXERCISES 


Chap. 1 

1. 10101101100111. 
2. 2(11)67. 

3, (10)(27)7. 


4. (a) The digits 2, (11), 6,7, expressed in the scale of 2—allowing 
for the (legitimate) addition of some 0’s—are respectively 
10, 1011, 0110, 0111. 


(b) The digits (10), (27), 7, similarly expressed are 1010, 11011, 
0111. Compare Ans. 1. 


5. 

123 456 456 -423- 
61 912 δ δ 
30 -1834- 114 -492- 
Ι5 3648 57΄ 984 
7 7296 28 -|968- 
3. 14592 14 -3936- 
Ι 29184 7 7872 
agp 3 15744 
56088 1 31488 
56088 


6. Since 707 = 7 x 101, apply the ‘s+ 1’ rule after separating 
N into pairs of digits, thus: 45,53, 31. Then 31 + 45 — 53 = 
23 = remainder. 


148 Answers to the Exercises 
Chap. 2 
1. (@) X 13 and add. (5) Χ 30 and subtract. 
21010101 21010101 
13 30 
2101023 2100980 
39 240 
210141 20769 
_13 270 
21027 1806 
91 10000000 180 
2193 9837131 7 
39 162869 69801 
258 3 7 
104 488607 = ἡ = 488607 
129 =3 x 43. 
2. 110999999 (997 = 1000 — 3) 
3 
11299999 
ie 
1329999 
3 
332999 
a 
33899 
9 


3989 Therefore the remainder is 1, 
9 and the quotient is 
998 111333 + 1 = 111334. 


997 
1 


3. Since 167 x 6 = 1002, use a similar procedure to the above 
but in this case add the multiples of 2. 


Answers to the Exercises 149 
Chap. 3 
1. 3/41 = -07317 29/41 = -70731 

7/41 = -17073 30/41 = -73170 

13/41 = -31707 


2. As with all recurrent cycles an algorithm—usually several—can 
always be found. I give the following because it differs markedly 
from those described in the text and shows something of the 
variety to be met with in this field. 

If m takes the successive values 1, 2, 3,...m, then 
11(2" — 2) defines the sequence 0, 22, 66, 154, 330, 682, 1386, 
etc. Dividing each term by 102" and summing, we have: 


ina & Seu 
002267573696145124... 


11(25 — 2) is, of course equivalent to 22(2"-!— 1); the 
latter part of this will become quite familiar. 


Chap. 4 

1. 3-1416. 

2. By definition F, = F3 — Fo 
Fo = Fa — Fo 
Fs = Fs — Fa 


Fa-1 = Fa+i — Fa 
Fn = Fa+e2 — ἔρωτι 


150 Answers to the Exercises 


Adding all these equations together we have on the left hand 
side the sum of the first Fibonacci numbers, whilst the terms on 
the right sum to ἔμ: — Fe. 


Since Fe = 1, the required sum is therefore ἔμ: — 1. 


Chap. 6 
1. 1001001001 x 111 = Jie 
100010001 x 1111 = Jie 


The factors common to both numbers are therefore those 
of Jie omitting those of 15 and /4, (3, 11, 37, 101). Hence 
they are 7, 13, and 9901. 


2. (See Tables 8 and 9) 
The sequence 91, 9091, 909091, 90909091, etc. occurs in the 
factorisations of J¢, Jio, 114,. .. Jax+1. Since 90909091 
does not appear among the prime factors of Jig it must there- 
fore be composite. 


3. 712 = 3333 X 33336667 = 3.7.11.13.37.101.9901 
Since 3333 = 3.11. 101, 
then 33336667 = 7.13.37.9901. 


Chap. 7 

1. The eight numbers relatively prime to 20 are 1, 3, 7, 9, 11, 13, 17, 
and 19 and therefore by Euler’s generalisation n can take any 
of the values 20k + 1, 3, 7, 9, etc. 
Thus, for instance 78 — 1 = 5764800 
and 198 — 1 = 16983563040. 


2. The index 12 appears at the following places in the table: 
m = 13, 21, 26, 28, 36, 42. Collectively these numbers contain 
the primes 2, 3,7, and 13 and it follows that these are all 
factors of 512 — 1. 

(5 being relatively prime to 12). 


3. 402 = 2.3.67. 
(402) = 402(1 — 1/21 — 1/3)(1 — 1/67) 
= 402 x 1/2 x 2/3 x 66/67 
= 6 x 1/3 x 66 =2 x 66 = 132. 


Answers to the Exercises 151 
Chap. 9 
1, Sx = 1001 — 9y 


Remainders on dividing by 5, 1 4 
Multiplying the integers 1 2 3 4 

by 4, removing multiples of 5;4321 

1 is in the fourth position, hence y = 4. 
Then 5x = 1001 —36=965. x= 193. 


2. For z= 26 = μι + πϑ, the only positive integer values of m 
and m are 5 and 1. 
Then x = m? — "ἢ = 24 
and y = 2mn = 10 
i.e. 242 + 102 = 26°, 


Chap. 10 

1. (a) 1. (6) 5. (Extrapolated from Tables 14a, δ). 

2. Divide each term as it appears, by 7; the remainder multiplied 
by three then provides the next term. In this way we get the 
sequence 1, 3, 2, 6,4, 5, 1,3, etc. The remainder thus appears 


as 1 when n=O, 6, 12,...= 6k and therefore 3% — 1 
is always divisible by 7. 


3. 5 and 8. (See Table 145.) 


Chap. 11 
1. 99999999 = 108 — 1 
Now 103 = —1 (mod 101) 
Therefore 108 = 1 
or 108 — 1 = 0 (mod 101) 
2. 25 = 32 = 9 (mod 23) 
210 = 81 = 12 
2u=24= 1 


222 = 1. or 222 — 1 = 0 (mod 23). 


3. Dividing throughout by 7 we find that 2908456 has the re- 
mainder 5, which is not a quadratic residue of 7. (δ) is 
therefore not a square. Dividing a,c, d,e, by 11, (6) has the 
remainder 2, and hence is not a square. 

Ans. (δ) and (e). 


152 Answers to the Exercises 
Chap. 12 
1. 29. 16733)35699 
33466 
2233)16733 
15631 
~ 1102)2233 
2204 
29)1102 
87 
232 
232 
“Ὁ 


2. It is easily seen that 3, 7, 11, 13, are not divisors so we set up 
ΠΡ) = 17.19.23.29 = 215441. 
17<p<29 
215441)1181027 
1077205 


103822)215441 
207644 


7797)103822 
7797 


25852 
23391 


2461)7797 
7383 


414)2461 
2484 


Therefore 23 is a factor. — 23 


Answers to the Exercises 153 


3, 10 4 65 
N = 1093709 
100 


204 937 
816 


2086 12109 
12516 
~ 407 
2093 = (2 x 1046) + 1 
2500 = 502 
Therefore N = 10472 — 50? = 1097 x 997. 


Chap. 13 

1, 15871 = 59 x 269. 
15853 = 83 x 191. 
15863 = 29 x 547. 


2. 13333 = 67 x 199. 
548497 = 53 x 79 x 131. 
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